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Abstract of 
Ultrasonic Absorption Along the Principal Axes of 
Single Crystalline Benzene 
by 
Alfred Erwin Victor 


Ph.D., Brown University, June, 1972 


Ultrasonic absorption measurements were performed as a function of 
temperature and frequency along the principal axes of single crystalline 
benzene. Measurements were made with longitudinal sound waves between 
6 and 26 MHz on all three axes and up to 34 MHz and 38 MHz on the a and 
c axes, respectively. The temperature range of the experiment extended 
from 170°K to 250°K; the range was limited by the characteristics of the 
transducer bond. 

The values of sound absorption (a) obtained (% 0.1 oo cm) 
are much higher than normally observed in solids (~ 10" higher than 
quartz). The magnitude of a corresponds to that obtained in previous 
experimental investigations of organic solids; however, the frequency and 
temperature dependence of a depart significantly from that which has 
been previously observed and from that predicted by theory. Present 
Eneories predict a to be proportional to the square of the frequency 
(v) in the temperature and frequency range of this experiment; however, a 


2 
was observed to depart from the v dependence as the temperature was low- 
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ered, the change being different for each axis. Relaxation effects were 
observed for each axis and indicate relaxation times of approximately 
lo”? sec; this is much slower than theoretically predicted (% aoe = sec). 
Specifically, for the a-axis the relaxation peak is centered at 17 MHz 

at 240°K and shifts to 12 MHz at 170°K. For'the b-axis, the relaxation 
Beak occurs at 26 MHz at 240°K and shifts to 20 MHz at 170°K. The c- 
axis data indicates the presence of multiple relaxation effects. It 
appears that one relaxation mechanism is occuring in the vicinity of 

14 MHz and another around 120 MHz. 

By comparing the change in the observed relaxation strengths as a 
function of temperature it was qualitatively determined that the sound 
absorption results from an exchange of energy between the lattice vibra- 
tions and the low frequency molecular oscillations. This was done by 
calculating the specific heat contributions to the relaxation strength 
for each internal mode for the upper and lower temperatures of the 
experiment and comparing the change with the observed values. Exam- 
ination of the molecular arrangement within the unit cell of benzene 
indicates that one should expect sound energy, which is propagated 
along each principal axis, to interact with the inter- and intramolecu- 


mar vibrations in a different manner. 
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Jigs INTRODUCTION 


The inherent complexity of molecular substances has been the stimulus 
for a large number of theoretical and experimental investigations. In 
many cases molecular structure has been precisely determined through die 
use of X-rays, neutron scattering, and Raman and infrared spectroscopy. 
The spectroscopic experiments together with the predictions of group theory 
have led to the determination of the vibrational spectrum of molecules, 
which, in turn, has been used to determine the force constants between 
neighboring atom-atom bonds. As the phase of a molecular substance changes, 
the structure of the molecule itself essentially remains intact. This has 
been borne out by Raman and infrared spectra experiments where the width 
of the lines may be altered but the fundamental intramolecular vibration 
frequencies remain about the same.” Additional low-frequency vibrations 
appear in the solid state which originate from the vibrations between the 
molecules in the lattice, i.e., the intermolecular vibrations. Since the 
intermolecular vibration frequencies are low, the coupling between the 
molecules in the lattice must be much weaker than the intramolecular couplings. 
For example, in solid benzene at0°C the highest lattice vibrational fre- 
quency is 3.15 x 10!2 Hz É whereas the lowest intramolecular frequency is 


1.215 x 1013 Hz.” 


In the harmonic approximation of lattice vibration theory 
the two types of vibrations are assumed to be independent of each other; 
however, when one takes into account the third and higher order terms of 

the intermolecular potential energy expansion, interactions between the 


two types of vibrations can take place. Raman and infrared spectra give 


experimental proof that such interactions do occur and are manifested in 
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combination frequencies which are observed but not predicted by the selection 
rules based on the harmonic Ee 

The possibility that the interaction between the intra- and intermolecu- 
lar vibrations could cause a large absorption of sound energy was first pro- 
posed by eterna: He hypothesized that the absorption would be a result 
of the slow exchange of energy between the two modes if their energy differ- 
ence was not too great. Using a theory which is based on a model of the 
molecular interaction in the gaseous state in combination with the Debye 
lattice continuum theory, he shows that (1) the resonance absorption is very 
large compared to other absorption mechanisms in solids, (2) that it is pro- 
portional to the square of the frequency, and (3) it is almost independent of 
explicit temperature variations. Later fie tenants substantiated his theore- 
tical prediction with absorption measurements on single crystalline benzene. 
The absorption values obtained were anomalously high (104 higher than crystal- 
line quartz) and could not be attributed to the classical absorption mechan- 
isms in solids; however, the values were of the same order of magnitude as 


EE elaborated on 


predicted by the resonance absorption theory. Rasmussen 
Liebermann's theory and made measurements on single crystalline cyclohexane 
and carbon tetrachloride. His theory and experiments qualitatively corrobo- 
rated Liebermann's results. The details of Liebermann's theory and Rasmussen's 
modifications of it will be discussed in Section 2.1. 

Further experimentation in the field was performed by Yun and Beyer” ?7° 
on single crystals of naphthalene and p-dichlorobenzene. The absorption 
values for both substances exhibited a frequency-squared dependence and 
little variation with temperature as ostensibly predicted by Liebermann's 


theory. Although the absorption values were relatively high for the fre- 


quency range investigated (5-95 MHz), they were about one order of magnitude 
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lower than those for benzene and cyclohexane. As a result, Yun and Beyer 
examined the contributions of the thermoelastic effect and dislocation 
damping more closely. It was found that the former mechanism was two 
orders of magnitude smaller than the observed absorption values. To deter- 
mine the contribution due to dislocation damping, a uniaxial pressure was 
applied to the crystals in small increments in a direction perpendicular 

to the direction of sound propagation. Although such treatment should 
increase the number of dislocations, no apparent change in absorption was 
detected. 

Wilson and Yun} extended the measurements in naphthalene down to 
liquid nitrogen temperatures. (Previous measurements in naphthalene were 
made between 200°K and 300°K.) At approximately 170°K a relaxation effect 
begins to appear. The absorption increases as the temperature is decreased 
below 170°K and peaks at a lower temperature; the position of the peak 
depends on the frequency. As the temperature is further decreased to the 
lower experimental limit (77°K), the absorption begins to decrease once 
again. Wilson and Yun interpret their results in terms of Liebermann's 
theory which does predict a relaxation effect at lower temperatures. How- 
ever, the predicted relaxation depends very strongly on the thermodynamic 


properties of the substance under consideration; consequently, the quantita- 


tive agreement between experiment and theory in this case may be fortuitous. 
The same calculation for benzene does not predict as strong a relaxation ef- 


fect although this work does indicate the existence of relaxation effects in 
benzene as well. 
: 12.13 . 
Danielmeyer extended the frequency range of the absorption mea- 


surements up. to the 100-MHz region. He made measurements in liquid and 
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solid bromine, benzene, and cyclohexane; however, his work had two short- 
comings. First, the thermal contraction of the experimental apparatus 
limited the temperature to -20°C; and second, the measurements on solids 
used polycrystalline samples, which made it difficult to correlate with 
previous single crystalline data. Danielmeyer did attempt to relate his 
measurements to previous investigations by proposing a theory which departs 
significantly from Liebermann's. The theory is based on a multiphonon 
interaction process and will be discussed in Section 2.2. 

All the absorption measurements discussed above were made on specimens 
whose crystalline orientations were unknown; consequently, the effects of 
mode conversion between the longitudinal and transverse waves which are 
present under such conditions practically rules out correlation with the 
actual structure of the substance under consideration. Heseltine, Elliot, 
and Wilson?! performed the first ultrasonic measurements on an oriented 
molecular crystal when they measured the longitudinal and shear sound 
velocities of single crystalline benzene. Since benzene is the basic 
aromatic hydrocarbon, many other experimental investigations have been 
performed on this substance; consequently, most of its structural and 
thermodynamic properties are known. It seemed appropriate, therefore, to 
carry out an experimental investigation on the ultrasonic absorption along 
the principal axes of crystalline benzene as a function of temperature. 


This is the subject of the present thesis. 





2. THEORY 


2.1 Liebermann's Theory 


E erman?” 


proposed that the chief cause of the large acoustic 
absorption in single crystals of organic materials is a relaxation process 
that converts, with a time delay, acoustic vibrations of the lattice into 
oscillations involving the internal degrees of freedom. That is, the transfer 


of energy from acoustic branches of the lattice spectrum into the optical 


branches is contemplated. He assumed that the absorption a could be expressed 


in terms of the standard equation for a relaxation process in a se 
a = a Sete te (2.1) 
ás o BE 1 + w*t? 


where s is the sound velocity, C; the internal specific heat per mole 
associated with the internal degrees of freedom, Se and O, are the specific 
heats per mole at constant pressure and constant volume, respectively. Tt is 
the relaxation time of the relaxation; its reciprocal tl is an approximation 
to the initial time rate of energy transfer between the external and internal 
degrees of freedom. Hence, a theoretical calculation of 1 yields a(w) if 
the other parameters in Eq. 2.1 are known. 


: 17 
Van der Waals forces are the main binding forces in a molecular crystal. 


They can be represented in various ways; the Lennard-Jones potential is 


a een (ans) 12 
sey = one» E 


where e and r_ are the Lennard-Jones parameters for, respectively, the 


convenient: 


depth of the potential well and the lattice distance from the reference 
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molecule,where U(r) is zero. Liebermann applied Eq. 2.2 to nearest. 
neighbor pairs of rigid or almost rigid molecules. Close to its minimum 
value, U(r) is nearly parabolic and therefore to first approximation it 
implies simple-harmonic oscillations, whether classical or quantum in 
detailed description. Ge of U(r) from a parabolic form can be 
treated as quantum mechanical perturbations that distort the zero-th order 
quantum mechanical wave functions. In particular, an acoustic signal 
applied at t = 0 can excite the lowest linear oscillator state of the 
lattice. Then, conventional use of first-order time-dependent perturbation 
theory will predict the time rate at which higher quantum states of lattice 
oscillations are excited. 

The Lennard-Jones potential deviates in a skewed fashion from parabolic 
form. Therefore, to a good approximation, only excited states of odd parity 
(in zero-th order) will be excited from the ground state, which is of even 
parity. When Liebermann applied first-order perturbation theory to cal- 
culate the transition rate for benzene, the third excited state was the 
highest energy eigenvalue contained in the resulting expression. He chose 
this state to be approximately equal to the lowest energy state of the 
optical branch. (Spectroscopic a available at the time would have 
shown this assumption to be invalid.) The calculation of the acoustic 
absorption coefficient is ostensibly completed by use of SE 2.1 under the 
assumption that the relaxation is from the third excited state of a longi- 
tudinal acoustic branch to the ground state of the lowest optical branch. 
When Liebermann calculated the value of a for benzene he used values of 
C; given by Lord, Ahlberg and Andrews. For substances for which the 


value of C; are not available in the literature, the Einstein model 





has been used to obtain the specific heat” 


SR Ee O3) 


R is the molar gas constant and x = hv/k,T where h is Planck's constant, 
v the frequency of the optical branch in question, Ko the Boltzmann con- 
stant, and T the absolute temperature. Rasmussen >> and Yun and Beyer” 
employed Eq. 2.3 to calculate the values of C; for substances used in 

their experiments. 

There are several assumptions, implicit or explicit, within the analysis 
for benzene. (a) The energy eigenvalue of the third excited state in the 
acoustical branch and an energy eigenvalue of an optical branch are essen- 
tially egual. (b) The center of mass of each molecule executes simple 
harmonic motion about its equilibrium position. The amplitude Xo of 
the associated lattice vibration of the acoustical branch caused by an 
applied acoustical signal is given by the classical equipartition 


of energy 


MX "wi = kel (2.4) 


where M is the molecular mass and w is the angular frequency of the 
lattice vibrations; i.e., the energy of the applied acoustic wave is 

already in thermal equilibrium with the lattice vibrations. (c) The time 
rate of transition in the relaxation process (171 in Eq. 2.1) equals the 
time rate for excitation of the third excited state in the acoustical branch. 
(d) The density of states in the acoustical branch is given by the Debye 
theory of specific heats of solids. 

] 


Liebermann's final prediction for 17 is 





ANERE T)” 
E O B 
T = Be; 2.5) 
16h(41*M) ™51%01"p 
where No is Avogadro's number, mol is the apparent internal vibrating 
mass, V is the atomic vibrational frequency associated with the transi- 


01 


Exon Of an optical branch from its ground state to its first excited state, 
VD is the Debye frequency, and A is the coupling parameter which is 
derived from the first anharmonic term in the Lennard-Jones potential. 

Liebermann” made measurements in single crystalline benzene at 6.4 MHz 
and 10.7 MHz at a temperature of 273°K. He found a proportional to wi 
for these two values and a nearly independent of temperature down to 
215°K. Upon substituting his values of a in Eq. 2.1 along with his experi- 
mentally measured sound velocity and the known values for the specific 
es, Liebermann found qt = 0.5 x 10”? sec. Substitution of known values 
for benzene in Eq. 2.5 yielded 1t = 1x 107? sec. This agreement tends to 
substantiate assumption (c) above. Also, either value of 1 makes 
wW T? << 1 + in agreement with the w*-dependence of the measured abscrption 
coefficients. 

Treating acoustic absorption in single crystalline cyclohexane (which 
is face centered cubic, fcc ), Kassa followed Liebermann's approach 
but attempted to improve on the details in three main ways. First, he 
replaced Liebermann's Lennard-Jones potential between nearest neighbor 
molecules by a similar potential between pairs of atoms in the nearest 
neighbor molecules. In doing this he assumed that each atom site had the 
mass Of one carbon and two hydrogen atoms and calculated effective 


Lennard-Jones parameters between the atom-atom pairs. Second, he assumed 





no orientation effect, i.e., the plane of the molecule was considered to 
lie parallel to the x-y plane of the crystal unit cell. Third, he replaced 
the Debye lattice spectrum by that proper for an fcc crystal (Brillouin*>). 
Fourth, he replaced the classical expression in Eq. 2.4 by the analogous 


quantum mechanical expression 


202 
=> 1 + A (2.6) 


272 
12k, T 


hw 
wcoth er 


| 2k 
B 


With these changes Rasmussen's expression for the relaxation time is 
2y 4 3 
o os. eee ee e E (2.7) 
01 2 In Ol m 


7 2 3 2 
Jm D oi Nu (un<M) 12(k T) 


where D is the coupling parameter analogous to A in Eq. 2.5 and pa 
is the maximum allowable frequency in the first Brillouin zone. The first 
Bern in Eq. 2.7 is nearly identical to Eq. 2.5 except that Vp has been 
replaced by Van and the numerical constants are different. The second 
term in the brackets in Eq. 2.7 turns out to contribute a small correction 
factor (approximately 5%) to the overall value of rl; consequently, 
though the basic assumptions between Liebermann's and Rasmussen's theories 
are different, the functional forms of the final results are essentially the 
same. Rasmussen's predicted values of absorption for cyclohexane were about 
50% higher than his observed values. His experimental measurements ranged 
from 3.8 to 8 MHz and showed frequency-squared dependence. Also, little 
change in absorption was observed with decreasing temperature down to 180°K. 
Considering the number and the vagueness of Liebermann's approximations, 
the agreement between his theoretically predicted and experimentally ob- 


served absorption values is fairly good. However, as Liebermann himself 
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points out, to describe the absorption process more accurately a more realis- 
tic intermolecular potential must be used to calculate the relaxation time. 
Although the Lennard-Jones potential was probably the best approximation 
available at the time (1959), it has been mainly used to describe the gaseous 
and liquid states; therefore, the Lennard-Jones parameters do not account for 
the interaction of molecules within the rigid structure of a solid where the 
relative orientation and intermolecular distances of the molecules are fixed 
compared to those in the fluid states. Rasmussen attempted to correct for the 
latter discrepancy by treating the atom-atom interactions between neighboring 
molecules; however, he assumed that the adjusted Lennard-Jones parameters 
which he derived for the nearest neighbor intermolecular atom-atom pairs had 
the same value whether the two atoms were at adjacent ends of two molecules 
or at greater separations. 

More recently D. E. Williams? has empirically derived an intermolecular 
potential for nine different molecular crystals (including benzene) which 
takes into account the actual atom-atom distances between neighboring mole- 
cules as determined from X-ray and neutron scattering experiments. Hence, the 
orientation of the molecules within the crystal structure is implicitly taken 


into account. The potential has the form 
E = Aee Ä Soe D + Bal explora) 
t Aog É Ton + Boy 2 exp(-Copt oy) 
+ Ay ) GE t+ Buy ) exp(=C ri) (2.8) 


and r represent the interatomic distances between the 


where Der: Toy HH 


carbon-carbon, carbon-hydrogen, and hydrogen-hydrogen atom pairs of nearest 
neighbor molecules. The maximun distance used in the summation is 6A°. 


The A, B, and C's are adjustable constants which are determined for each 
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type of atom-atom interaction by correlating the potential energy with 
known values of the thermodynamic and structural properties for the sub- 
stance under consideration. Bernstein”? used Eq. 2.8 with Williams! values 
of the A, B, C constants to calculate the values of the acoustical and 
optical phonon frequencies for crystalline benzene. The agreement between 
the calculated and experimental values is goad, which suggests that Eq. 2.8 
would be useful in future theoretical investigations of acoustic absorption 
in molecular crystals, provided that the crystalline structure is known in 
detail. 

Liebermann's calculation for the transition probability implicitly 
assumes that the number of transitions per unit time from the third excited 
state of the acoustical branch to the ground state of the lowest optical 
branch is equal to the number of transitions per unit time for the reverse 
process. Such a supposition is questionable from statistical mechanical 
considerations. The Bose-Einstein distribution predicts the occupation 
number of the lower energy state to be greater than that of the upper state 
which implies that the number of transitions per unit time should be unequal 
for the two processes. 
si 


is proportional to v A A 


Equation 2.5 indicates that Tt D D 


is considered to be a constant, then +t would decrease with increasing 
temperature. Also, the values of Cs en directly proportional to the 
temperature whereas the sound velocity is inversely proportional. There- 
fore, in the region where wt << 1, i.e., where a is proportional to w*, 
the temperature dependences of the parameters in Eq. 2.1 counteract each 


other. Thus, a is predicted to increase only slightly with increasing 


memperature.. A. O. Williams” points out, however, that Yo is directly 
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proportional to the sound velocity through the relation?” 


Y = con}? (2 9) 
D unV avg : 


(A similar expression exists for ye Or theme icc ee Subse- 
quent to Liebermann's theoretical formulation, velocity measurements as a 
function of temperature have been made in eyalohexane, - carbon tetrachlo- 
m 8 9 E 9 13 f 

ride, naphthalene p-dichlorobenzene, and benzene, and in all cases the 
sound velocity exhibits a strong linear dependence on T; s decreases as 
T increases. (In benzene the velocity changes by 14% between 170°K and 
Zo0°K. ) Bee showed that when the experimental velocity values are 


-l from Eq. 2.5, a is predicted to fall very fast with 


used to obtain T 
increasing temperature in the temperature range of the experimental inves- 
tigations (180°K-310°K). This is inherently implausible and in strong 
disagreement with the approximate independence of a on T from the mea- 
surements of Liebermann, Rasmussen??? and Yun and Revenue 

Finally, reference to experimental spectroscopic data on crystalline 
benzene upsets the energy eigenvalues used by Liebermann [assumption (a) 
above]. At 273°K the highest energy level in the acoustic branch is found 
Ro correspond to 105 em”! (O = 150°K)* and the lowest in the optical branch 
to 406 cm”? (0 = 5820K) 18 | a ratio considerably exceeding three. Conse- 
quently, higher energy states of the acoustical branch ee considered 


in order to conserve energy between the two levels. This approach will be 


Ronsidered in the next section. 


os 


2.2 Danielmeyer's Theory 


AS proposed a theory which would be applicable to the 
acoustic absorption in any molecular crystal, providing the Raman and 
infrared spectra, the specific heats, and the elastic constants of the sub- 
stance are known. He postulates that the sound absorption is primarily a 
result of multi-phonon processes between the acoustical and optical energy 
levels, and for some substances between different levels of optical branches. 
In using the multi-phonon concept, the transition probability is calculated 
between the different energy levels involved which results in a relaxation 
time for the process. The relaxation time is related to the acoustic ab- 
sorption by deriving an expression for the complex sound velocity by com- 
bining (1) the equation of motion for a solid, (2) the equations of state 
for a solid as described by Liebfried and GE and (3) the reaction 
rate equation of Herzfeld and Se The resulting equation for the com- 
plex sound velocity is then separated into its real and imaginary parts 
which, in turn, yield expressions for the velocity dispersion and sound 
absorption, respectively. Both quantities are described in terms of a 
relaxation strength relating the ce and thermodynamic properties of 
a substance and the relaxation time, which combines the transition proba- 
bility between energy levels with the anharmonic properties of the inter- 
molecular potential. To derive the above expressions, Danielmeyer assumes 
the Debye-Einstein model for the molecular solid which stipulates that the 
frequency of each optical branch is constant throughout the first Brillouin 
zone and that the frequency of each acoustical branch is linearly propor- 


tional to the appropriate sound velocity (longitudinal or transverse) 
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throughout the first Brillouin zone, with a maximum frequency of vp at the 
outer edge of the zone. Furthermore, only transitions between the longitudi- 
nal acoustic branches and the optical branches are considered because in most 
solids the longitudinal acoustic energy levels lie closer to the optical 
levels; hence, a greater transition probability would be expected. 

To conserve energy in a phonon-phonon interaction, the sum of the ener- 
gies of the phonons which are destroyed must equal the energy of the phonon 
which is created. Since the optical levels are usually much higher than the 
acoustical levels, Umklapp processes (U-processes) play an important role in 


such a transition. Using the Debye-Einstein model, Ziman shows that the 


: ; : S T 
fraction of U-processes involved in a three-phonon interaction is 
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where ve is the frequency of the optical branch, Vp is the Debye frequency, 
q is the magnitude of an arbitrary acoustic wave vector, and s is the sound 
velocity. Equation 2.10 shows that the number of U-processes increases as 

the difference between ” and Vp increases. ER and vp are proportional 


to the Einstein temperature Op and the Debye temperature On» respectively. 
Hence, Eq. 2.10 implies that the ratio of U-processes increases up to the 
limit where 0,./9,, = 2, at which point no two acoustic wave vectors would be 
large enough to satisfy the energy condition and still be less in magnitude 
than the maximum allowable wave vector in the first Brillouin zone, dt 

In general the ratio 0,,/0, is greater than two for most organic solids; 
consequently, it is necessary to consider higher-order phonon processes in 
the calculation of the transition probability. In benzene, for example, 


0/0, = 3.85 ; therefore, at least four acoustic phonons must combine 
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to produce one optical phonon; hence, a five-quantum process is required. 

To calculate rigorously the transition probability of multi-phonon 
processes (m-quantum processes) one must use the higher order, anharmonic 
terms in the potential EE Bach anharmonic term used yields an addi- 
tional contribution to the transition probability. However, a generalized 
intermolecular potential for organic solids was not known at the time; con- 
sequently, Danielmeyer proposed that the m-quantum process could be con- 
structed from a series of three-quantum processes. The transition proba- 
bility for each three-quantum process is calculated using first-order 
time-dependent perturbation theory and the first anharmonic term in the 
potential energy expansion for an isotropic solid. The energy levels 
between the initial and final states are considered to be virtual energy 
states which are short lived, and consequently notwell defined as a result 
of the uncertainty principle. Danielmeyer shows that for a given m-quantum 
process there will be (m-2) three-phonon processes involved. 

To determine the order of the phonon process involved for a particular 
substance one must know the energy relationships between the various acous- 
tical and optical levels. This information can be readily obtained from 
Raman and infrared spectroscopy data. In some cases the separation between 
energy levels in the optical spectrum is large enough to consider two or 
more separate optical groups. This occurs in cyclohexane, for example, 
where the spectroscopic data show one group of optical levels falling in 
the range Op = 333°K to 751°K and a second group between Op = 1154°K to 
2094%K, The difference between the closest energy levels of the two groups 


is 402°K whéreas the difference between the highest acoustical energy level 
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(characterized by On , which is 130°K for this case) and the lowest optical 
level is 203°K. Therefore, two relaxation processes are possible. The 
acoustic-optic process with 0./0) = 2.55 requires the annihilation of 
three acoustic phonons to create one optical phonon, a four-phonon process. 
The optic-optic transition with 97/09) = 3.1 requires four acoustical 
phonons plus one optical phonon from the lower optical group to be anni- 
hilated to create one optical phonon in the upper optical group, a six- 
phonon process. In general, the optic-optic transition will influence the 
sound absorption only when it is of a higher order than that of the acoustic- 
optic transition because a larger value of qt is predicted for the higher 
order process and, as will be shown later, in the region where wt << 1, 

a is proportional to wir, Hence, the optic-optic transition will be taken 
into consideration only when the process is at least one order higher than 
the acoustic-optic process for the same material. 

The levels within an optical group are spaced so close together in 
comparison to the energy spacing between groups that, once energy is coupled 
into the group, the excess energy is rapidly distributed among the levels 
within the group. In many cases where a vibrational relaxation has been 
observed in a molecular gas or liquid, it has been found that the total 
vibrational specific heat can be associated with a single relaxation time; 
consequently, the entire optical group can be represented by a single char- 
acteristic internal temperature, E Danielmeyer assumes that this 
generalization holds true in the molecular solid and uses it in deriving 
an expression for the acoustic absorption (the derivation will be discussed 


below). 


Danielmeyer examined the spectroscopic data of all the materials on 
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which acoustic absorption measurements had been made, i.e., benzene ,º 

7 9 E 9 e SS 
cyclohexane, naphthalene, p-dichlorobenzene and bromine, and found the 
p-dichlorobenzene was the only other substance besides cyclohexane for 
which the optic-optic transition was important. (Note: All equations 
which are not directly referenced below have been taken from Danielmeyer's 
original work’? ) To treat the spectroscopic, data in a general manner, the 
ratio 97/07 is defined in terms of a transition parameter 


Bro o Vara 
AS 2 J == l J (2:11) 
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where j= j' + 1 and the notation Op is replaced by SR The index j 
denotes the optical level under consideration and is equal to 1, 2, or 3 
corresponding to the bottom of the first optical group, the top of the 
first optical group, and the bottom of the second optical group, respec- 
tively. The transition parameter is greater than two for all the organic 
substances discussed; consequently, an m-quantum process must be considered 
for each case. 

The quantitative calculation of the m-phonon transition probability 
incorporates a number of simplifying assumptions. (1) The m-phonon process 
can be described as the sum of a series of three-phonon processes. ‘There 
are (m-2) three-phonon processes for a given m. (2) Each acoustical and 
optical energy level is considered to SSC, simple harmonic oscillator. The 
isothermal relaxation time between oscillator states <f| and |i> is 
given by 


eb = - -12 
Tr Peg Pee (2 ) 


where Tr is the isothermal relaxation time. Pig is the probability that 


a particle in energy level |i> will make a transition to energy level 
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BR). Pe; is the corresponding probability for the reverse reaction. 


(3) The number of particles in an energy state j is given by the Bose- 


Einstein distribution 


1 
n; = oA > (2:13) 
ei e J EEN 1 


For acoustic phonons the high temperature approximation of Eq. 2.13 is 


assumed, 1.e., T >> 0 hence, the occupation number becomes 


D 5 
p SH 

Nai ~ hv (2.14) 
ai 





where the index ai refers to the ith acoustic phonon. 
(4) The number of particles making transitions via a three-phonon 


process from state |i> to <f| is given by the product DË. e which is 


related to the standard transition probability expression derived from 


first-order time-dependent perturbation theory by 


E a 2 5 
Wc = Hc or | LN pla) att, SEKR (218) 


Here H. is the matrix element <£|H,li> and H, is the perturbation 


fi 


Hamiltonian obtained from the first anharmonic term in the intermolecular 
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potential. The density of states is assumed to be given by the Debye 


approximation and is expressed in the form 
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where V is the volume and Ji is the ith acoustic phonon wave vector. 
The delta function ensures energy conservation between the two states in 
question. By combining Eqs. 2.11 and 2.16 and using the properties of the 
delta function to integrate over the volume, Eq. 2.15 can RR E in 


terms of the transition parameter A: 





me 


= Blue o 2 
Caner ce qn He. | Ca) 


where So is the sound velocity in the frequency range in which no dis- 
persion exists. 

(5) The een Hamiltonian is obtained from the theory of vibra- 
Bons in an elastic continuum. In the classical theory of elasticity it 
is assumed that the potential energy density at any point is a quadratic 
function of the local strain, i.e., Hooke's Law. The first anharmonic 
potential energy density term involves a cubic function of the local strain 


and is given by the expression” 


1 
Ha = EN ) SE S ES dv TTE) 
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where the A.. are rrelated to the third order elastic con- 
13k%mn 
stants and the S; F are the displacement gradients. Equation 2.18 
3 


is usually written in a more convenient form utilizing the 
creation and annihilation operators. After several substitutions involving 
the relationships between the displacement gradients, normal coordinates, 
and the creation and annihilation operators, Eq. 2.18 can be rewritten for 


E : 32 
the special case of a single three-phonon process as 
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Here the =" 5nd a are the creation and annihilation operators, 
SE SÉ 


respectively. The delta function ensures momentum conservation within the 
first Brillouin zone (Normal processes) and allows for the possibility of 


E . e 
U-processes through the reciprocal lattice vector g. B is defined as the 
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coupling parameter and is related to a summation of the E For 
longitudinal sound waves propagating in an isotropic solid B can be 


expressed in terms of the Grtineisen constant,” Y: 


B 
Er (2220) 
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m is the second-order elastic constant associated with the longitudinal 


sound velocity. The Grüneisen constant, in turn, can be related to the 





macroscopic thermodynamic properties of a solid bei" 
rs alle : (2-21) 
pC x 


where ß is the volume expansion coefficient and e is the isothermal 
compressibility. Thus, the coupling parameter can be related to the bulk 
properties of the solid. Substitution of Eq. 2.19 into Eq. 2.17 allows one 
to calculate Weg and Wes for the process being considered, i.e., 
optic-optic or acoustic-optic. The transition between two optical levels 


> > 
involves one acoustical phonon (a, = 9,12 and two optical phonons 


> > > 
(4,7 = do» duer P ds). By use of Eqsememls, 2.17, and 2.22, Eq. 2.12 


becomes 
mk Tv 3B2 
= Do, = nz aya [2 [6/7 a]? = intro, 
d 242p9s 3 
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where the quantities K and J, are defined as indicated. (J, will be 


O arises from consideration of the 


discussed below). The factor (e 
occupation numbers of the phonons at equilibrium as given by the Bose- 
Einstein -distribution. (Note that Liebermann neglected to take this into 


account, Which was one of the basic criticisms of his theory). 


The same equations are used to calculate the three-phonon acoustic-optic 
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process, with one exception. Because the process involves two acoustical 
phonons, the delta function in Eq. 2.22 can be satisfied by many values of 
Ee and Eë consequently, the integral over the density of states must 
be modified. Looking ahead to the m-quantum process in which T acoustic 
phonons are involved, Danielmeyer defines a density-of-state function Jp 
in terms of the transition parameter A. -He then calculates Jp for several 


values of T and from the set of values obtained he empirically derives a 


function J,. J, turns out to have the simple form 


A A 
J, = Do for ı>1ı 
(2. 24) 
J = À? for A<I 
PT and A are related by the expression 
pas 5 Reg (2.25) 


where in the argument of the cosine function 0 and 1 are used for the 
acoustic-optic and optic-optic transitions, respectively. The transition 


parameter A, then, is the fundamental quantity of the calculation. 
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For a three-phonon acoustic-optic process e becomes 
-1 = - 
5 on 10 
un2y Dette! "e "TM. E 
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where the index 0 refers to the top of the acoustical band. Once again 
the exponential factors arise from consideration of the occupation numbers 
in the two states. Equation 2.26 does not fully agree with Eq. 2.23, but 
een" has attempted to keep their main structures alike. 
Assumption (1) above states that for an m-quantum process there are 


(m-2) contributions to the total transition probability. Also, for a given 
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substance it was pointed out that the optic-optic transition would influence 
the absorption only if the process was at least one order higher than the 
acoustic-optic transition. With these criteria, Danielmeyer shows that the 
number of contributions to e is related to the number of acoustic 
phonons, IT , involved in a given m-quantum process. For the optic-optic 
transition 
(m-2) = P (2,274) 

whereas for the acoustic-optic transition 

(m-2) = T-1 (2.27b) 
Equations 2.27 indicate that there are T-1 and T-2 intermediate virtual 


energy states for the respective transition types. To calculate E 


for 
the m-quantum process, the transition probability from the initial real 
energy state to the first virtual energy state must be calculated; then the 
transition probabilities between the successively higher virtual energy 
states are calculated and finally, the transition probability between the 
highest virtual state and the final real energy state is calculated. When 


this is done, the total transition probability for the optic-optic m-quantum 


process is 
2 
= ES en (2.28) 


and for the acoustic-optic case it is 


2 
oe = wl? mT-1)T (,0/T_,-9/T), (2.28) 
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Comparison of Eqs. 2.28 and 2.29, respectively, with Eqs. 2.23 and 2.26 


show that they differ only via the factors of [ involved. The factor K 


is raised to a power equal to the number of three-phonon processes that 
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contribute to the total relaxation time, i.e., T or T-1, and the entire 
equation is multiplied by the same factor. The density-of-states integral 
Jr is evaluated by using Eqs. 2.24 and 2.25. 


To relate the relaxation times to the sound absorption, Danielmeyer 
derives an expression for the complex sound velocity in an isotropic 


molecular solid. The exchange of energy between the internal and lattice 


vibrational modes is described by the well known relaxation rate equation”? 
(T; - T,) 2.30) 


T, , the characteristic temperature of the internal vibrations, is assumed 
to represent all the internal vibrations within the optical group under 
consideration (see discussion on page 16). T, is the temperature that 
accompanies the sound wave due to the localized compressions and rarefac- 
tions; therefore, it is periodic in nature. In addition to the localized 
temperature variations, the sound wave will induce localized deformations 


in the solid which are described by the equation of motion” 


Oe pee (2.31) 
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where NEG: is the derivative of the stress tensor SE with respect 
to the jth coordinate; 5; is the second time derivative of the particle 
displacement in the ith direction; p is the density. Under the assumption 
that the deformations are small enough to allow neglect of second-order 


: : 37 
effects, s, is related to the strain tensor RE by 


l1 
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The perturbation of the thermal and mechanical properties of a solid that 


is caused by the passage of a sound wave can be described in terms of the 
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thermal and caloric equations of state (Liebfreid and Ludwig”) which relate, 


respectively, the pairs of quantities T, and Ir and Ze and Së 
to variations of the Helmholtz free energy of the system, which, in turn, 


can be expressed in terms of macroscopic thermodynamic quantities. For an 


isotropic solid the thermal equation of state becomes 
Tee 
T YC, p(T_-T) 
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where aaa is the isothermal elastic constant associated with any longi- 
tudinal sound velocity and y is the Griineisen constant which can be cal- 
culated from Eq. 2.21. The time derivative of the caloric equation of 


Beate is 


i d a d ° 8 
Cc. SE (T,-T) + C; a + TYC Es; = 0 (2,34) 
where cm and a are the specific heats associated with the internal 
(optical) modes and acoustical modes, respectively. Assuming that the 
i(ker-wt) . 
sound wave propagates as a plane wave described by e with ampli- 
tude small enough so that T. and T can be represented by small de- 
partures from the equilibrium temperature, Danielmeyer combines Eqs. 2.30 


through 2.34 to obtain an expression for the complex sound velocity 
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The real part of Eq. 2.35, which describes the velocity dispersion yields 


z — pa (2.36) 
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The imaginary part, expressed in terms of absorption per wavelength, is 
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The quantity r is the relaxation strength, defined by 
(C_-C_)z ci 
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E. in Eq. 2.38 is related to the mechanical properties of the solid by 
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where E and Se are the adiabatic elastic constants for an isotropic 


solid. The isothermal relaxation time Tr has been replaced by E Berck 
relaxation time at constant pressure and entropy. The two are related by 


the expression 


1 A (2.40) 


Danielmeyer uses Leg instead of Tp SO that the absorption values obtained 


in the solid can be compared with those of the gaseous and liquid state 
where Es is the relaxation time customarily used to describe the relax- 
ation process. °° 

Equations 2.35 and 2.36 are in the standard form of the relaxation 
theory for liquids and gases; however, they have been obtained from a 
solid state viewpoint rather than from the fluid state approach used by 
Liebermann. The properties of the solid are explicitly contained in the 
relaxation strength, i.e. Es > and the relaxation time, i.e., the coupling 


parameter B. 
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In the region where NE << 1 , Eq. 2.37 can be written in terms of 


a as 
2n2r1 Eve 
a, = — == (2.41) 
O 
For wt >> 1 
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Equation 2.41 indicates that the absorption coefficient is proportional to 


the square of the frequency as in Liebermann's case. Also, a, is directly 


dependent on r and E Tee is inversely proportional to temperature, 


whereas r is directly proportional; therefore, the two temperature-dependent 
effects tend to counteract each other and result in a slight temperature 


dependence of a, , the absorption decreasing with decreasing temperature. 


x 


Equation 2.42, however, indicates that a is independent of frequency. 


i 


Also, since both r(l-r)7! and KS are directly proportional to the tem- 


perature, a, is influenced by temperature changes a great deal more than 
a 

Using values available at the time for parameters contained in the 
equations leading up to and including Eq. 2.41, Danielmeyer calculated the 
absorption coefficient for single crystalline benzene at 6 and 10 MHz at 
270°K to be 0.099 cm”! and 0.310 cm, eis: These values compared 
quite well with Liebermann's experimental values of 0.09 cm”! and 0.24 cm”! 
which were measured at 6.4 and 10.7 MHZ, respectively, at 273°K. However, 
Danielmeyer's theory does not agree as well with other experimental results. 
A comparison will be made below. 


From Eq. 2.37 , the frequency at which the sound absorption 


coefficient (per wavelength) reaches a maximum is 
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Equation 2.43 predicts the maximum absorption per wavelength in solid 
benzene at 275°K to occur at 800 MHz. In liquid benzene the corresponding 
frequency has been theoretically predicted as approximately 600 MHz, Y 
Experimental values of Heasell and Lamb" and Hunter and Dardy > qualita- 
tively corroborate the predicted value although the measurements have been 
made only up to 510 MHz. Hence, Danielmeyer concludes that since no relax- 
ation effects appear in the liquid state until the frequency reaches about 
600 MHz, solid benzene will exhibit a similar behavior and relaxations will 
not appear until the frequency reaches the vicinity of his predicted 
800 MHz. It will be shown in this work that although the frequency-squared 
dependence of the absorption in solid benzene is approximately true down to 
250°K, the experimental absorption values exhibit relaxation effects at 
much lower frequencies than Danielmeyer predicts; the amount of departure 
from the v* dependence is a function of crystal orientation. 

As in Liebermann's case, Danielmeyer's relaxation times (Eqs. 2.28 and 
2.29) are strongly dependent on the Debye frequency contained in the factors 
Kl and (T = 4.52 for benzene). All of Danielmeyer's predicted values 


assume v, independent of temperature. If 1, is calculated for benzene, 


D T 


using Eq. 2.29 for the temperature range of 170-250°K with Vy constant, 


Te increases with decreasing temperature and the absorption in the region 
where wt << 1 is predicted to be practically independent of temperature. 


Conversely, if v, is assumed to vary in accordance with Eq. 2.9 (using 


D 


published velocities of BE Tr decreases by one order of magni- 


tude in the same temperature range which, in turn, predicts a corresponding 
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decrease in a. Although this does not agree with published experimental 
results, the temperature dependence of the absorption values obtained in 

this work does show a decrease of absorption with decreasing temperature, 
at rates depending on the crystal orientation and the ultrasonic frequency. 


Therefore, in contrast to the effect of a temperature dependent Vy on 
Liebermann's predicted relaxation times, Danielmeyer's equations still 
predict a plausible result. The exponential factors contained in 
Eqs. 2.28 and 2.29 and the velocity term, Se OS ofiset the 
dependence of Vp with temperature. 

Table I shows a comparison of the experimental and theoretical 
relaxation times for both Liebermann's and Danielmeyer's theories. The 
experimental values were obtained by substituting the measured absorption 


values into Eq. 2.1 and 2.41, respectively. Equation 2.40 was used to 


convert the Ven obtained from Eq. 2.41 into Tape The theoretical values 
were calculated by the various investigators as indicated. When Daniel- 


) for naphthalene and p-dichlorobenzene, he used 


meyer calculated (Te ae 


the absorption values given in the original paper of Yun and Beyer; - 
however, they later published an Erratum?” in which they noted that these 
experimental results as originally reported were in error due to a miscal- 
culation. Danielmeyer apparently did not see the erratum; consequently, 
he was misled into believing that his theory predicted an arden of mag- 
nitude agreement with experiment for these two substances. Danielmeyer's 
calculated values for sde in Table I are shown in parentheses for 
naphthalene and p-dichlorobenzene. If these values were correct, Daniel- 


meyer's theory would agree favorably with experimental results except for 


cyclohexane. The values of Tere calculated from the corrected 
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absorption values are shown without parentheses. It can be seen that 
agreement is rather poor for all substances except benzene. Conversely, 
Liebermann's theory indicates an order-of-magnitude agreement with all 
the experimental results to date; however, as pointed out earlier some 
of the fundamental concepts of Liebermann's theory are questionable. 
Although Danielmeyer attempted to account for some of the weaknesses in 
Liebermann's theory, the experimental results do not corroborate his 
predictions. 

Unquestionably, the theoretical analysis of acoustic absorption in 
molecular solids is an extremely complex problem. Danielmeyer and 
Liebermann attempted to simplify the problem to a point where a quanti- 
tative calculation could be made; however, the results of both theories 


indicate that further theoretical development of this subject is desirable. 
2.3. Absorption Due to Thermal Conduction 


The experimental absorption values are the sum of the vibrational 
absorption and the "classical" absorption due to thermal conduction. The 
latter quantity can be calculated from the theory of Woodruff and 
Ehrenreich. ‘> In the region where wt << 1 , the absorption due to thermal 
phonons is given by 

q E Cie (2.44) 
ps“ 
where «x is the thermal conductivity and y is the Grüneisen constant., 
Instead of using experimental values of «x , Danielmeyer substitutes the 


e. H 
expression 
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in Eq. 2.44 where d is the lattice constant. Equation 2.44% becomes 


2 
a = Zu" (2.46) 
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By using Eq. 2.46 to calculate ang > it can be shown that in none of 


the cases considered does D exceed 1% of the experimental absorption 


values. 
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3. Experimental Apparatus and Procedures 


Set Introduction 


The properties of organic crystals differ considerably from those of 
ordinary "hard" solids; consequently, the handling of organic materials poses 
a number of onerous problems for the experimenter. Benzene freezes at 5.2°C 
and has a high vapor pressure of 36 mm Hg at that temperature which causes 
the solid to sublime very rapidly unless it is kept in a sealed container. 

To provide a work area that can be maintained below 5.2°C a refrigerated 
glove box which could be kept between -10°C and -15°C was constructed. pene 
zene crystals are relatively soft; therefore, special care had to be taken 
during cutting and facing operations, to avoid straining them. Benzene vapor 
is extremely re consequently, precautions had to be taken to exhaust 
the working area and to avoid breathing the fumes generated during prepara- 
tion of the crystal for measurement. Hence, the glove box was located in a 
hood room and an organic vapor respirator was worn. 

Since high values of acoustic attenuation were anticipated, double 
transducer operation would have been desirable; however, because of the rapid 
sublimation rate, a special transducer bonding technique had to be used which 
sealed the sample from the atmosphere while at the same time bonding the 
transducer onto the sample. Double transducer operation was incompatible 
with the bonding technique. 

Prior to perforining ultrasonic absorption measurements on oriented 
benzene samples a number of auxiliary experimental techniques had to be 
developed and/or applied. Liquid benzene was purified by vacuum sublimation. 
Large single crystals of benzene were grown from the melt. Orientation was 


accomplished by use of a polarizing microscope and a universal stage. Cutting 
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and polishing of the oriented sample was accomplished with a specially 
designed mechanical apparatus. Finally, the sample faces were sealed (the 
transducer bond provided the seal on one face) and ultrasonic measurements 
were made by the standard pulse-echo technique. Measurements were made from 
6 to 26 MHz on all axes; additional measurements up to 38 MHz were made on 
the a and c axes. The characteristics of the bond restricted the tem- 
perature range between 170°K and 250°K. Since measurements were made as a 
function of temperature, only relative absorption values were obtained. 
Attempts were made to ascertain the attenuation due to the transducer bond; 
however, only an approximate estimate was obtained. Each step of the experi- 
mental procedure will be briefly described in the following sections; the 


specific details are given in the appendices. 


3.2 Vacuum Sublimation Technique 


The success of single crystal growth is largely dependent on the purity 
of the liquid when it is initially put into the crystal-growing apparatus. 
The benzene used in this experiment was purchased from the Fisher Scientific 


% Pure'' and free of thiophene, the 


Company and was certified to be ''99 Mol 
principal contaminant usually found in benzene. However, it was still neces- 
sary to degas the liquid benzene of dissolved air prior to crystal growth. 
Among the many methods available for distillation, the vacuum sublima- 
tion technique ° is particularly adaptable to materials that have a high 
vapor pressure in the solid state. Also, the procedure is carried out in a 
completely closed system; consequently, the problem of toxic vapors is elimi- 


nated. A vacuum sublimation system, constructed as an integral part of the 


crystal-growing apparatus, was designed to allow relatively large quantities 
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of benzene to be degassed fairly quickly. 

The vacuum sublimation process is based on the fact that a substance 
maintained at a pressure or temperature below its triple point values cannot 
exist as a liquid. Phase transitions occur from the vapor to the solid, or 
vice versa, depending on whether heat is liberated or absorbed. On a phase 
diagram of pressure versus temperature (p-T diagram), the solid-vapor equil- 
ibrium curve (vapor-pressure curve) depicts the point at which a phase tran- 
meron will take place for a given set of conditions. The substance will be 
a solid if both pressure and temperature are below the triple point; however, 
if the pressure is kept constant below the triple-point value and the tem- 
perature is raised, the substance will sublimate when the temperature reaches 
a value at which the solid and vapor are at equilibrium, i.e., when the tem- 
perature coordinate on the p-T diagram intersects the vapor pressure curve. 

During a sublimation process, molecules leave the surface of the solid 
(the sublimand) and either remain in the vapor phase or recombine with the 
sublimand. Hence, the sublimation rate is directly proportional to the sur- 
face area of the sublimand and inversely proportional to the amount of recon- 
bination that takes place. The surface area is controlled by the physical 
dimensions of the vessel which contains the sublimand; therefore, the vessel 
is made as large as is practicable. Recombination can be prevented by pro- 
viding an external means to carry off the molecules in the vapor before they 
can recombine. The external means in this case is a vacuum which is main- 
tained throughout the system. 

To dispose of impurities in it, the vapor must be recondensed to a solid 
at a temperature and pressure at which the impurity melee ties remain in the 


vapor phase while the molecules of the substance to be purified are condensed. 
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If the pressure of the system is constant for both the sublimation and 
condensation processes, and is below the triple points of all the substances 
in the vapor, the condensation temperature must be above the triple point 
temperature of the impurities but below the triple point value of the 
substance to be purified. 

The triple points of oxygen and nitrogen are approximately 2 mm Hg/54°K 
and 96 mm Hg/63°K, respectively. 17 The triple point of benzene is 
36 mm Hg/2780x 38 Ineneronpe, ir frozen benzene which contains dissolved air 
is sublimed at a pressure below 2mm Hg, the oxygen and nitrogen will pass 
off as gases. If the ensuing vapor is subsequently condensed at 77°K (at a 
pressure below 2mm Hg), the oxygen and nitrogen will remain in the vapor 
phase and continue to pass through the system, whereas the benzene will re- 
main in the solid phase, freed of oxygen and nitrogen impurities. 

In the construction of the vacuum sublimation system, it is necessary 
to insure that the components used in the system do not contaminate the sub- 
stance which is being processed. Since benzene is a solvent for most greases, 
any stopcocks or glass joints which require vacuum grease as a lubricant 
cannot be used. Benzene also attacks most rubber materials used to make 
O-rings; however, viton is an exception and can be used safely. With these 
criteria, the system was fabricated from glass tubing of various diameters. 
Teflon stopcocks and glass-to-glass O-ring joints (with viton O-rings) were 
used where required. 

A schematic drawing of the system is shown in Fig. 1. The teflon stop- 
cocks are designated by E the glass-to-glass O-ring joints by Ji and 
the normal vacuum system valves are designated by Ve The supply tube, the 


storage tube, and the large tube of the sublimate trap were made of glass 
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tubing with an outside diameter (0.D.) of 52 mm. Since the sublimation rate 
depends a great deal on the surface area of the sublimand, these tubes were 
made as wide as possible yet narrow enough to fit into a standard one-liter 
dewar flask. The glass manifold and the inner tube of the sublimate trap 
were made from 10-mm O.D. tubing. The length and diameter of the inner tube 
in the sublimate trap determine the capacity of sublimate that can be deposi- 
ted during the process. The wider the tube, the less likely it is to become 
logged. The inner tube length is two-thirds that of the outer tube. Thus, 
approximately 100 ml of sublimate could be deposited at a time. The actual 
procedure used to sublimate the required amount of benzene and the method 
for putting the degassed benzene into the crystal growing apparatus are 


described in Appendix A. 


Bes Crystal Growth 


During the orientation procedure (Section 3.4) and the subsequent 
preparation of the benzene single crystals for ultrasonic measurements 
(Section 3.5), a significant part of the original piece of crystal was 
trimmed away because of the cutting and melting operations involved. It 
was therefore necessary initially to grow large pieces of single crystalline 
benzene. The final oriented sample was 1.9 cm in diameter and 1.5 cm long. 
The crystal-growing apparatus and procedure used were adopted from a method 
developed by Liebermann and Beer The crystals are grown from the 
melt in a sealed container, by using a cold copper probe as the nucleation 
point for the crystal seed. Once the seed is formed, the crystal is grown 
by slowly varying the temperature of the melt, the probe, and/or the ambient 
air. The actual temperatures and their rates of change vary for each type 


of material in ways that must be determined principally through extensive 
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trial-and-error experimentation. Once the procedure is established, large 
pieces of single crystal can be grown (typically 10 x 10 x 5 em). 

The fundamental condition for crystallization to occur from a melt is 
that the total free energy of the solid-liquid system be decreased. For 
each volume of crystal that is created there is a decrease in the specific 
free energy of the solid and an increase in its surface free energy, the 
latter being a function of the surface tension at the solid-liquid interface. 
Mer Chere is a critical size for the nucleus and a concomitant critical 
change in free energy which must be established for the volumetric effect to 
outweigh the surface effect. When this condition is fulfilled, the nucleus 
will support stable crystalline growth. SE shows that for a given 
substance the critical change in free energy is decreased by (1) supercooling 
the liquid and (2) providing an external surface for nucleation. 

Paes” shows that, in general, the rate of crystal growth is inversely 
proportional to the shear viscosity of the liquid and directly proportional 
to the amount of supercooling. He also points out that the temperature 
distribution within the solid-liquid system is a major factor in the quality 
of crystals grown and is one that can readily be controlled experimentally. 
If one considers the liquid-solid interface as a plane surface, the heat flow 
equation per unit area is 


dT dT 


eo se L 
K Ta = AH yf + SA GES (3.1) 


Here Ke and K, are the thermal conductivities of the solid and liquid and 


dT /dx and dT, /dx the corresponding temperature gradients. AH is the 
change in the specific enthalpy, i.e., the latent heat of fusion per unit 
Volume, and f is the rate of linear growth in the x direction. Equa- 


tion 3.1 shows that as the crystals grow in the x direction the latent heat 
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of fusion must be dissipated; this requires that the solid-liquid interface 
be cooled as the crystal grows. The heat can be dissipated through the solid, 
provided the proper temperature gradients exist. Ideally, the temperature 
gradient of the system should be such that a minimum temperature is at the 
center of the solid mass, a positive gradient exists between the center of 
the solid and the solid-liquid interface, and a positive gradient (smaller 
than that in the solid) exists from the interface to the outer extremity of 
the system. If the heat of fusion is not dissipated completely, the inter- 
face may become slightly warmer than some part of the liquid farther out in 
the melt. The result is a negative temperature gradient in that part of the 
melt, which leads to growth of unusable crystals. 

In general, the intrinsic properties of surface tension, viscosity, and 
thermal conductivity of a liquid determine the difficulty of crystal growth 
for a given experimental apparatus. Experimentally, one has control over 
the temperature gradients and their time rates of change in the system and, 
to some extent, the degree of supercooling of the liquid. The presence of a 
cold probe in the liquid enhances nucleation of the seed crystal and also 
provides a path for the dissipation of the latent heat of fusion. 

The entire crystal growing apparatus is built into a commercial refrig- 
erator and is shown in Fig. 2. A large copper block (3/4 in, thick by 7.5 in, 
long by 6 in. wide) was made to attach onto the cooling coils of the refrig- 
erator. Attached to the block is a tapered, stainless steel lid assembly 
which is used to seal a small pyrex bell jar. SE probe is threaded 
through the copper block, passes through the lid, and extends into the jar. 
Two glass tubes and a thermocouple probe also pass through the lid-and-block 


assembly into the jar. A plexiglass box containing several heater elements 
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Fig. 2b. Crystal growing jar and associated equipment. 
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Fig. 2. Crystal growing apparatus. 
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encloses the entire apparatus and is used to control the ambient air tempera- 
ture that surrounds the jar. The temperatures of the probe, the surface of 
the lid, the ambient air, and the liquid in the jar (the melt) are controlled 
my electric heaters to accomplish crystal growth. 

The details of the copper block and stainless steel lid assembly are 
shown in Fig. 3. The large thermal mass of the block minimized temperature 
fluctuations in the copper cooling tip. The lid is tapered to ensure that 
the cooling tip is the only nucleation point in the liquid and its surface 
temperature is always kept above the freezing point of benzene by an internal 
heater element. The latter condition precludes crystallization of benzene 
apor on the lid which, if allowed to occur, would interfere with the 
initial stages of crystal growth on the probe. 

The copper cooling probe is threaded through the copper block and passes 
through a teflon bushing inside the lid (see Fig. 3). A viton O-ring in the 
Bushing seals the area between the probe and the bushing from the atmos- 
phere. A copper-constantan thermocouple and an electric heating element are 
embedded inside the probe. The part of the probe going through the teflon 
bushing has a 1/2-inch diameter and the tip of the probe was machined to a 
hemisphere of radius 1/4-in. The probe tip was polished to a fairly smooth 
surface with emery cloth and Alumina polishing compound. 

One of the glass tubes which passes through the block-lid assembly is 
Msed to fill the jar. The filling tube goes through the freezer compartment, 
through the top of the refrigerator, and terminates outside the refrigerator 


Meethe joint J, shown in Fig. 1 (see also Fig. 2a). The other glass tube 


T 


Which passes through the block-lid assembly passes from tne freezer compart- 


ment through the right side of the refrigerator and is connected to the 
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vacuum system. The area between the outside surface of each tube and their 
respective passage in the lid~block assembly is sealed by a quick-disconnect 
O-ring seal (CAJON Ultratorr Fitting) which is soldered into the top of 
the copper block. 

The thermocouple probe which passes through the lid-block assembly is a 
1/8-inch 0.D. stainless steel sheath which has an iron-constantan thermo- 
couple embedded inside. The probe is six inches long and can be positioned 
within the melt or inside the lid at any desired vertical level. The outer 
surface of the probe sheath is also sealed with a quick disconnect O-ring 
connection. 

The crystal growing jar (see Fig. 2b) is a small pyrex bell jar with a 
Capacity of approximately 600 ml (purchased from the American Scientific 
Company, catalog number JR-4875). The jar is set in a heater element which 
conforms to the contour of the jar base. The heater element is in turn set 
on a plywood base 3/4-in. thick. Four threaded rods which are attached to 
the copper block pass through the wooden base, and the jar is sealed to the 
lid by tightening the wing nuts on the rods at the bottom of the wooden base. 
All seals in the jar-lid assembly are made with viton O-rings and without 
any vacuum grease. The jar heater is controlled by a thermistor sensing 
element which is taped on the outside wall of the jar and is thermally 
insulated from the ambient air. 

The air temperature around the jar assembly is controlled by a grid of 
heating elements which are threaded through the lower sides of the plexiglass 
box (see Fig. 2b). A thermistor extends down from the wooden base and is 
positioned in the middle of the heater grid; its control unit will be 


described below. A linearly polarized light source, which is used to check 
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visually for flaws during crystal growth, is attached to the inside back wall 
of the box; when the box is in place, the light source is directly behind the 
center of the jar. The four threaded rods pass through the bottom of the box 
and allow securing the box in position. A copper-constantan thermocouple, 
inserted through the center of the right side of the box, is used to monitor 
the air temperature at that point. 

There are four heaters in the apparatus that require control. Power to 
the probe and lid heaters is provided by Variacs to establish the tempera- 
ture gradient within the cooling probe. The air and melt temperatures and 
the temperature gradients in the liquid are controlled by half-wave SCR (sil- 
icon control rectifier) proportional controllers in conjunction with the 
thermistor sensing elements. The SCR units were designed and built in this 
laboratory by Cuomo”- and Bliniek after the original design of scot 

The temperature gradient in the liquid had a radial component due mainly 
to the ambient air temperature, and a vertical component due to the melt 
heater. It was found that once all the heater controls were set to their 
optimum position (to be described later), only the melt-heater control had 
to be varied during crystal growth. A motor drive unit, attached on the melt 
temperature SCR controller, continuously decreased the melt temperature at a 
very slow rate. The fastest rate compatible with good crystal growth was 
Found to be 0.12°C/hr; it was achieved by turning the SCR heater control at 
trate of 0.36 rpd. 

All temperatures were monitored with a Mosley 7100B strip-chart recorder 
which could be read to within +0.005 mV on its most sensitive scale. Because 
of the large thermal inertia of the copper block, the probe and lid tempera- 


ture did not exhibit temperature fluctuations in excess of +0.1°C when at 
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equilibrium. Fluctuations in the melt temperature and the air temperature 
were approximately +0.25°C and +1°C, respectively. 

Initial experiments with the apparatus showed the need to incorporate 
slight modifications in the equipment described. It was found that the best . 
crystal growth was obtained when the tip of the copper probe penetrated 
approximately one inch into the liquid. In this position, the lid tempera- 
ture did not much influence the probe temperautre. However, so much surface 
area was in the liquid that three or four single crystal seeds would usually 
form and this limited the ultimate size of a piece of single crystal that 
could be grown. Consequently, a teflon collar was put around the probe tip 
to reduce the area available for nucleation. Subsequently, one or, at most, 
two single crystal seeds would cone 

An additional jar heater was wrapped around the top of the jar (see 
Fig. 2b) and was controlled by the same power source as that for the lid 
heater. RS served to prevent crystallization on the jar walls immediately 
below the outer extremity of the lid; the lid heater was not able to do this. 
Considerable improvement in crystal growth was realized thereby. 

It was originally intended to use the stainless steel thermocouple probe 
to monitor the melt temperature as it was changed; however, it was found 
that when the melt temperature reached approximately 3.5°C, crystals started 


to form on the probe which interfered with the main crystal on the copper 


probe. Therefore, once the amount of supercooling (=1.5°C) and the optimum 
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It is believed that a cone shaped probe with a small tip radius would 


overcome the problem to some extent for the tip of the cone would present 


. 7,54 
an even smaller nucleation area. See Rasmussen. 
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Fig. 4. Single crystal of benzene immediately 
after being removed from crystal growing jar. 
(Approximate dimensions 7 x 7 x 10cm). 
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rate of decrease in the melt temperature (0.12°C/hr) were determined, the 
thermocouple probe tip was positioned inside the lid to monitor the lid 
memperature. 

Since all the temperatures are interdependent, it was necessary to 
determine experimentally which combination of Variac and SCR controller 
settings yielded the best crystal growth. The actual procedure that was used 
and some problems associated with the apparatus are discussed in Appendix B. 
The method used to cut the crystal and prepare it for optical orientation is 
described in Appendix C. Fourteen crystals were grown with the apparatus; 
each one yielded crystals large enough for two or three samples. A piece of 
Single crystal, which had just been removed from the jar, is shown in Fig. 4. 


Its approximate dimensions were 7 x 7 x 10 cm. 


3.4 Orientation of Single Crystalline Benzene 


Since it was intended to measure a in known crystallographic directions, 
the next step was to determine these in the crystal sample whose growth was 
described in Section 3.3. Benzene is an orthorhombic crystal; it therefore 
has three axes, "a", "b", and "c", which are mutually perpendicular but each 
with different length of unit cell. Attempts to determine the axial direc- 
tions were first made with the Laue back-reflection X-ray technique; however, 
these were unsuccessful, the interpretation of the complex X-ray patterns 
being excessively difficult. Fortunately, benzene crystal is transparent, 
colorless, and strongly birefringent (biaxial). Hence, methods of optical 
orientation with polarized light were available er, 

The orthorhombic structure happens to be the only one for biaxial 


crystals in which all three principal crystallographic axes coincide, respec- 


tively, with the three "indicatrix axes." In general, the biaxial indicatrix 
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is a triaxial ellipsoid in which the three mutually perpendicular semiaxes 
are equal in magnitude to the three principal indices of refraction, 

EK: ny < n,. Two optic axes and the acute bisectrix (the bisector of the 
angle between the optic axes) lie in the no n, plane, the optic plane. 
For the orthorhombic case, the acute bisectrix coincides with one of the 
axes of the crystal, and the optic plane coincides with one of the (100) 
planes. Since all three of the crystal axes are mutually perpendicular, one 
must lie in the optic plane, perpendicular to the axis defined by the acute 
bisectrix, while the third axis must be perpendicular to the optic plane. 
Mullens” has determined that for benzene the acute bisectrix and the optic 
plane correspond, respectively, to the b axis and the a-b plane of the 
eeystal. 

A beam of plane-polarized (p-p) light entering a birefringent crystal is 
split into two components that travel at different speeds. One of these 
components is perpendicular and the other parallel to the principal section; 
i.e., to the plane containing one of the optic axes and the propagation vec- 
tor of the light passing through the crystal. If the birefringent crystal is 
set between two polars of a polarizing microscope (i.e., a polarizer and an 
analyzer) that are maintained crossed at 90°, maximum extinction of the light 
will occur when the crystal is rotated to a position for which the planes of 
vibration of the light components in the crystal coincide with the planes of 
vibration of the polars. At this point a relationship between the microscope 
axes and the principal section has been established: one of the optic axes 
Will lie in a plane parallel to one or the other of the vibration planes of 
the polars. 


The polarizing microscope used to achieve this first result (as well as 
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the subsequent orientation operations) is a fairly conventional compound in- 
strument with the following auxiliary features: (1) The sample is held in a 
multiaxes stage which allows the crystal to be rotated about three mutually 
perpendicular axes. (2) A special objective lens (30X) is used which allows 
enough separation between the lens and the microscope stage to permit inser- 
tion of the multiaxes stage on the microscope stage. (3) There are perpendi- 
cular crosshairs in the eyepiece (ocular) which correspond to the directions 
of the vibration planes of the polars when the latter are held crossed at 
90°. One of the polars (the polarizer) is located below the microscope stage 
and the other (the analyzer) is located between the objective lens and the 
eyepiece. The two polars can be crossed at any arbitrary angle but are usu- 
ally used in the 45° or 90° position. (4) A removable lens (Bertrand lens) 
can be inserted between the analyzer and the ocular. The microscope can be 
used as a conoscope or an orthoscope depending on whether the Bertrand lens 
is inserted or removed, respectively.°° The conoscopic mode was used in this 
experiment. 

After the crystal is mounted on the multiaxes stage and placed on the 
microscope stage, a convergent beam of p-p white light is brought to focus 
on the bottom face of the sample. This beam diverges as it passes upward 
through the crystal; upon leaving, it EE through the objective lenses, the 
analyzer (crossed at 90° to the polarizer) and the Bertrand lens. The result 
is an interference fringe pattern which is focused at a point between the 
Bertrand lens and the ocular. The size and resolution of the pattern seen 
through the ocular depend on the magnifying power of the objective and the 


ocular. 
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Because the beam of polarized light diverges while passing through the 
crystal, some part of it will be parallel (or nearly so) to one of the optic 
axes. Along any optic axis, light of either polarization travels with the 
meme speed; hence, this particular direction will be shown by a black dot 
surrounded by fringes, because the polars are crossed. Superimposed are dark 
bands called isogyres which start from the optic axes as sharply defined 
(straight or curved) lines, but become progressively broader away from this 
axis. When the acute bisectrix of the crystal is aligned parallel to optic 
axis of the microscope, the isogyres form a cross which coincides with the 
crosshairs in the ocular. The two optic axes patterns are equally spaced 
about the center of the cross on one of the isogyres. When this pattern is 
observed with all the axes on the universal stage in their zero positions, 
Ene crystal is oriented with its b axis along the optic axis of the micro- 
Scope, its a axis along the isogyre containing the optic axes figures, and 
its c axis along the other isogyre of the cross. The acute bisectrix 


Brenn is shown in Fig. 5. In Fig. 5 the a axis is parallel to the E-W 





Fig. 5. Acute bisectrix pattern (from 
Hartshorne and Stuart [58]). 
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direction, the b axis is normal to the plane of the paper at the center of 
the cross, and the c axis is in the N-S direction. 

For a complete theory of the interference patterns observed with a 
conoscope see References 57 and 58. For a description of the manipulations 
involved in this experiment see Appendix D. The procedure used enabled 


Erientation to within +2°. 


3.5 Sample Preparation 


The methods used to prepare the sample for ultrasonic measurements None 
patterned after those of Heseltine and Elliot.” An oriented piece of crystal 
is cut into a cylindrical sample and frozen into an annular brass sample 
holder. The faces of the sample are lapped to make them flush with the flat 
sides of the sample holder and are sealed with thin pieces of mylar. The 
mylar seals prevent sublimation of the sample and on one face provides a 
bonding surface and a ground connection for the ultrasonic transducer. The 
sample holder is mounted onto a specially designed housing which incorporates 
the electrical connections necessary to excite the transducer and provides 
the means for positioning the sample holder in a dewar assembly which, in 
turn, enables one to make ultrasonic measurements as a function of temperature. 

All the sample preparation is done in the refrigerated cold box with the 
mechanical apparatus shown in Fig. 6. This apparatus is made entirely of 
brass and is used in a vertical or a horizontal position. The main components 
consist of the dovetailed track, the cutter-planer housing, the sample holder 
housing, the planer, two cylindrical cutting tools, and an L-shaped stand 


which supports the track in the vertical position. A breakdown of these 


Somponents is shown in Fig. 7. 





E 





Fig. 6a. Mechanical apparatus in horizontal position 
(Used to melt flat face on sample). 





Fig. 6b. Mechanical apparatus in vertical position 
(Used to cut cylindrical sample). 
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The cutter-planer housing and the sample holder housing fit onto the 
track and can slide along the track axis or be held stationary by tightening 
set screws in the base of each housing. Both housings are machined so that 
Eheir vertical edges are exactly perpendicular to the track axis; hence, 
parallel to each other. A piston which is two inches longer than the cutter- 
planer housing fits into a cylindrical hole in the center of the housing and 
is free to move along the track axis. A mounting plate on one end of the 
we ton enables either the planer or a cutting tool to be attached to the 
piston. The planer is a flat brass disc which is heated by a nichrome wire 
heating element embedded behind it. The cutting tools are thin-walled 
(1/64-in. thick) brass tubes (3/4-in,and 7/8-in, diameter) which are heated by 
a nichrome wire heating element wrapped around the upper end of each tube. 
Inside the sample holder housing there is a removable cavity (see Fig. 7). 
The cavity provides the means to attach the sample holder to the housing 
during the sample facing procedure and is removed during the sample cutting 
procedure. 

The sample holder (see Fig. 7) is a brass ring which has a 2-in,0.D., 
is 0.593 in. (1.54 em) thick, and has a 3/4-in. diameter hole bored through 
the center. Several additional holes are incorporated which are used (1) to 
mount the sample holder on the removable cavity, (2) to attach a cover 
plate over the mylar seal on the reflecting surface of the sample, (3) to 
mount the sample holder on the ultrasonic cavity (to be described later), 
and (4) to insert a thermocouple into the sample holder. 

The procedure used to cut and polish the sample with the mechanical 
apparatus is described in Appendix E. The final sample is 1.54 cm long and 


1.9 cm in diameter and has one face (the reflecting surface) sealed with a 





~ 


EE 


piece of 1/'%-mil thick mylar. The next step is to bond the transducer onto 


the opposite face of the transducer. 


3.6 Transducer Bonding Technique 


To obtain a good coupling of the sound energy from a transducer into a 
solid the bond between the two must be able to support the type of mechanical 
displacement induced by the transducer and be as thin as possible, the latter 
condition being necessary to ensure that a minimum of energy be lost in the 
bond itself. The bonding agents and procedures used generally depend ne 
nature of the bonding surface and the type of transducer being employed, 
ke., it is usually more difficult for a bond to support shear waves than 
longitudinal waves. 

The frequency dependence of the sound absorption yields useful 
information concerning the properties of the material being investigated. 

The upper frequency limit at which sound absorption measurements can be made 
is principally determined by the intrinsic attenuation of the sample, i.e., 
the higher the intrinsic absorption the lower the upper frequency limit. 

When measurements are made as a function of temperature, the relative values 
of a for different frequencies at the same temperature and the change in 
the relative values of a with changing temperature are of primary interest. 


When this type of measurement is being made, the most consistent results are 


usually obtained if a single transducer is used and is operated at its funda- 


mental frequency and at as many odd harmonic frequencies as is possible, 
given the intrinsic attenuation of the sample. Transducers with different 
fundamental frequencies used on the same sample can lead to erroneous results. 


(This point will be discussed later.) 





re, 
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Benzene is very soft and fragile; consequently, it cannot withstand a 
great deal of pressure before it becomes strained or cracked. Because the 
sound absorption in benzene is so high, a transducer with a relatively low 
fundamental frequency (2.0 MHz) had to be used to determine the frequency 
dependence of the absorption with temperature. Low fundamental-frequency 
transducers are inherently thicker than the higher fundamental-frequency 
transducers (5 MHz and above) normally used with solids and, in general, are 
more difficult to bond onto a sample. Hence, the fragile benzene surface 
and the need for a thick transducer necessitated the use of a peculiar bond- 
ing technique to obtain a usable ultrasonic signal. 

A brief description of the bond will be given here; the details are 
presented in Appendix F. The transducer is bonded onto a piece of 1/4-mil 


thick aluminum-coated mylar with Nonaq stopcock grease at room temperature. 


|The transducer-mylar combination is, in turn, bonded onto the benzene sample 


ı with a glycerine-water solution. The latter bond is made in the cold box at 


about -10°C. The mylar seals the sample from the atmosphere and provides an 

electrical path from the ground electrode of the transducer to the sample; 

the latter is at ground potential when the ultrasonic measurements are made. 
Once a bond was made, it would last for approximately five days before 

it began to deteriorate. It was found that a bond could be broken and 

remade with a fresh transducer-mylar combination two to three times. The 

low and high temperature limits of the ultrasonic measurements were restric- 

ted by the characteristics of the bond. The lower limit appears to depend 

on the linear expansion coefficient of the benzene along its three principal 

axes. Although it could not be precisely determined, it is believed that the 


bond broke as a result of the benzene-mylar combination pulling away from the 
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transducer due to thermal contraction. The linear expansion coefficient is 


greatest along the c axis? 


and the bond broke at a higher temperature along 
this axis than along the a or b axis. Consequently, measurements along 
the c axis were only taken down to 180°K whereas the lower limit along the 
a and c axes was 170°K. The high temperature limit was due to an apparent 
change in the bond characteristic when the sample temperature was raised from 
250°K to 260°K. Measurements were taken at intervals of 10°K, first with the 
temperature decreasing and then with the temperature increasing. If the tem- 
perature was allowed to go above 250°K after one cycle and subsequently 
brought below 250°K to initiate another cycle, it was found that the measure- 
ments of the two cycles were inconsistent. However, if an upper limit of 
250°K was used, the measurements of succeeding cycles were consistent to 
within 3-4%. 

An x-cut, coaxially plated quartz transducer with a fundamental frequency 
of 2MHz was used. The transducer has a ]/2-in. overall diameter and a 3/8-in, 
active diameter. Transducers with fundamental frequencies of 2.5 MHz and 
5.0 MHz were also employed; however, inconsistent results were obtained with 
the three transducers and are attributed to the difference in the transducer 
thicknesses. (This point is discussed more fully in Section 4.2 and 
Appendix G.) Since the transducer with the lowest fundamental frequency 


yields the most information within a given frequency interval, the 2 MHz 


transducer was used exclusively for the measurements. 


3.7 Ultrasonic Measurement Procedure 


The prepared sample is mounted onto the ultrasonic cavity as shown in 


Fig. 8. The ultrasonic cavity provides the electrical connections necessary 





-ajetd Zuryunou JeMop 
a4} OF payoejzze ÁJTABD DTUOSPAITON 


*6 wou] 


“(OT *3TI 998) 
ÁJTAPD DTUOSPAITON 






Bee 








259- 


to excite the transducer and is designed to fit into a dewar to enable meas- 
urements to be made as a function of temperature. The sample holder slides 
Over the four posts of the cavity base and rests on the bottom plate. Four 
coil springs, one on each post, provide spacing between the spring-loaded 
nen electrode and the transducer. The top plate is secured to the four 


posts and has a screw threaded through its center which allows pressure to 


foe exerted onto the electrode plate which, in turn, adjusts the pressure 


that the electrode exerts on the transducer. It was found that only a light 
pressure on the transducer was necessary to obtain a good ultrasonic signal. 
The ground connection to the transducer is provided by the aluminum coated 
mylar which is in contact with the brass sample holder, i.e., the whole 
cavity assembly except for the center electrode is at ground potential. A 
coaxial cable is connected to the electrode plate and there are access holes 
in each plate to allow a thermocouple to pass through and be inserted into a 
hole in the sample holder. 

The ultrasonic cavity is attached to the mounting plate of the dewar 
assembly with two stainless steel rods as shown in Fig. 9. The rods are 
first attached to the upper plate of the cavity and then screwed into the 
bottom of the mounting plate. The coaxial cable from the cavity is connected 
to the mounting plate and the thermocouple passes through an access hole in 
the top of the plate. 

A schematic of the dewar assembly is shown in Fig. 10. A liquid helium 
dewar was modified to allow the ultrasonic cavity to be inserted into a tem- 
perature controlled environment. A small brass can which is wrapped with 
heater tape is located in the lower part of the inner dewar. Two copper- 


constantan thermocouples are embedded between the can and the heater tape. 
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FIGURE 10. SCHEMATIC OF SAMPLE INSIDE OF DEWAR ASSEMBLY 








Pele 


One of the thermocouples is used as the sensing element for an SCR propor- 
tional controller which controls the current in the heater and the other 


measures the temperature of the outside of the can. The ultrasonic cavity 


is suspended in the can by the mounting plate which, in turn, is secured 


- 


to the upper flange of the dewar assembly. The outer dewar is filled with 


liquid N, which provides the refrigerant for the system. 


The temperature inside the can is controlled by a commercially manufac- 


'tured feedback control system. A block diagram of the main components of 


the system is shown in Fig. 11. (All components were manufactured by R. I. 
Controls, Inc., Minneapolis, Minn.) The Data Trak programmer unit supplies 
command signals to the Thermac Controller. The Thermac Controller provides 
the Power Controller with an input signal which is proportional to the dif- 


ference in the command signal of the programmer (corresponding to the desired 


‚system temperature) and the feedback signal from the system being controlled 


(corresponding to the actual system temperature). The Power Controller is a 
full wave SCR proportional controller which supplies the amount of current to 
the heater required to bring the system temperature back to the desired 
value. The Data Trak programmer gives either a steady state command signal 
for constant temperature control or a programmable, rate signal which changes 
the temperature as a function of time. With this system extremely stable 
temperature control and smooth rates of temperature change can be realized. 
When the temperature inside the dewar assembly is at equilibrium, temperature 
fluctuations would not exceed +0.1°K. The equipment is capable of even better 
temperature control; however, there is a considerable amount of thermal lag 
in the system because of the large amount of brass used in the mechanical 


apparatus. 
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Prior to transferring the sample from the cold box into the dewar, it 
is important that the temperature inside the can be stabilized at about -10°C. 
If the sample is subjected to a thermal shock, the bond may be destroyed 
and/or the sample may crack. (Note: The temperature at the thermocouple 
next to the surface of the can is normally lower than the temperature inside 
the can; consequently, a relation between the two must be determined experi- 
mentally.) Once the sample is in the dewar, the sample temperature is moni- 
tored with the copper-constantan thermocouple in the sample holder. The 
thermocouple is connected to a Mosely 7100B strip-chart recorder which can 
be read to within +0.005 mV (~ +0.1ºK). After the sample is initially put 
into the dewar, the temperature is allowed to stabilize for approximately 
15 minutes. When equilibrium is reached, the Data Trak programmer is turned 
on and the temperature is decreased at a rate of 5-8°K per hour. Upon reach- 
ing a desired temperature, the system is allowed to come to equilibrium for 
25-30 minutes prior to making any measurements. Measurements were taken 
every 10° between 170-250°K. (The lower limit along the c axis was 180°K 
AS indicated in Section 3.6.) 

A block diagram of the electronic system used to make the ultrasonic 
measurements is shown in Fig. 12. All the components except the oscilloscope 
(Fairchild, Model 766H) were manufactured by the MATEC Corporation of Warwick, 
R. I. The model numbers of the equipment used are shown in Fig. 11; in some 
cases more than one component is incorporated into the same unit. The master 
synchronizer triggers the oscilloscope, the pulsed oscillator, and the pulse 
comparator. The latter two are alternately triggered to allow the wave forms 
generated by each unit to be alternately displayed on the oscilloscope; the 


persistence of the cathode ray tube of the oscilloscope makes the waveforms 
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appear as if they are superimposed. The pulsed oscillator supplies an elec- 
trical pulse of short duration (% 5 usec) to the transducer which converts 
the electrical energy into acoustical energy with a frequency equal to the 
transducer fundamental frequency or an odd harmonic of the fundamental. 
The acoustical energy is reflected off the back face of the sample and is 
received by the transducer. The electrical signal generated by the trans- 
ducer is amplified and demodulated in the receiver and the echo pattern is 
displayed on the oscilloscope. 

The pulse comparator provides a calibrated comparator pulse which can 
be used to measure the difference in the attenuation between any two pulses 
in the echo train. The calibrated exponential generator provides an alter- 
nate method of attenuation measurement. An exponential waveform is super- 
imposed on the echo train and is adjusted to make the peaks of the pulses 
in the echo train fall on the exponential. This method provides a convenient 
means of checking the quality of a pattern (which should be exponential) 
which, in turn, can be related to the quality of the bond and the extent that 
diffraction and nonparallelism effects have on the attenuation. A more de- 
tailed description of the electronic components and their uses is found in 
the eene af? 

The specifications for the Hodel 950 low frequency plug-in, which is 
used with the Model 6000 Generator and Receiver unit, state the frequency 
range to be 1-20 MHz; consequently, for measurements above 20 MHz a Model 560 
Transmitter and Receiver was used (frequency range 10-320 MHz). In practice, 
it was found that acceptable patterns were obtained with the Model 950 up to 
| 26 MHz; consequently, it was possible to double-check measurements between 


10-26 MHz. Measurements with the two units usually agreed to within 2-3%. 
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The attenuation measurements were made as a function of frequency at 
each temperature. The measurements were made at a frequency of 6 MHz or 
higher, the upper limit depending on the orientation. No data were obtained 
at 2 MHz, for preliminary measurements showed the attenuation values obtained 
were inconsistent with those of higher frequencies. This is attributed to 
Ehe lange diffraction loss at this frequency which could not be accurately 
accounted for with the correction procedure used. (This point will be 
discussed further in Section 4.2.) 

For highly attenuating materials, there are relatively few pulses in 
the echo train; consequently, an accurate attenuation measurement with the 
exponential generator is not practicable. Therefore, the pulse comparator 
method was used in this case. A change in the relative value of attenuation 
with temperature is most accurately obtained by taking the difference in the 
attenuation between two echos in the pattern. The use of two pulses has the 
inherent advantage that any change in the electrical characteristics of the 
electronic components is automatically accounted for, i.e., a change in the 
electrical characteristics will have an equal effect on each pulse in the 
pattern and will not influence a difference measurement between two peaks. 

The mechanics of making a measurement with the pulse comparator are 
relatively simple. The height of the comparator pulse is matched to the 
height of the first pulse in the pattern. Then the comparator pulse is 
positioned next to another pulse farther out in the echo train and is atten- 
uated in calibrated steps until it is the same height as the second pulse 
selected. The amount of attenuation inserted can be read directly off the 
dials on the pulse comparator unit to within +0.1 db. Because the attenua- 


tion in benzene is so high, only two peaks were available for measurement in 
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the patterns for 14 MHz and higher; consequently only the first two peaks 
were used at all frequencies. Pictures of typical patterns obtained with a 
2 MHz transducer are shown in Fig. 13. The patterns shown depict the atten- 


uation along the b axis at 170°K and 250°K at 6, 10, and 14 MHz. 
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Fig. 13. Observed ultrasonic absorption patterns for benzene, b-axis, 
170°K and 250°K. (Scale: 2 volts/cm vertical and l0u sec/cm 
horizontal.) 
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4, Experimental Measurements 
4.1 Experimental Results 


Initially measurements were taken with transducers having fundamental 
frequencies of 2.0, 2.5, and 5.0 MHz. However, as mentioned in Section 3.5, 
the results obtained with the three transducers were inconsistent. Although 
the frequency dependence appeared to be the same for the three cases, the 
absorption values did not agree to within experimental error, i.e., the 
curves obtained on an a versus vV plot for each transducer would be dis- 
placed with respect to each other by a distance which fell outside the 
experimental error. This anomaly was first attributed to a bond loss; how- 
ever, further experimentation did not bear this out. Regardless of whether 
each transducer type was bonded onto different samples or successively 
bonded onto the same sample, the results were the same. Since so many dif- 
ferent bonds yielded the same discrepancy between the attenuation values, 
it was concluded that the bond characteristics were not the predominant 
factor causing the inconsistency. Although it was not proved conclusively, 
the discrepancy seems to be associated with the differences in the transducer 
thicknesses. Experimental measurements which tend to support this conjec- 
ture are presented in Appendix G. Because of this effect all the measure- 
ments presented here were taken with a 2.0 MHz transducer since it yields 
the most information within a given frequency interval. 

The data presented here were taken with five separate samples. One 
sample was oriented along the a-axis, two along the b-axis, and two along 


the c-axis. During each run the sample temperature was lowered in ten 


de 


During the course of the experiment measurements were taken on 21 samples. 
Various combinations of bonding techniques and transducers were tried before 
reproducible bonds could be made and good ultresonic patterns were obtained. 
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degree intervals to the minimum value. At each interval the temperature 

was allowed to come to equilibrium for 25-30 minutes; then the measurements 
were taken. After the minimum temperature was reached, measurements were 
taken in the same manner with increasing temperature up to 250%K, The cycle 
was repeated three times, which resulted in six sets of attenuation measure- 
ments at each temperature. 

The average sound absorption values for the a, b, and c axes, which 
have been corrected for diffraction effects, are tabulated in Tables II, 
III, and IV, respectively. Very weak signals were obtained at 30 MHz and 
higher; consequently, measurements could only be made to within +1 db 
(normal accuracy is +0.1 db). As a result, the data for 30 MHz and above 
are somewhat erratic and not as reliable as those for the lower frequencies. 
The standard deviation for each point is not greater than 7%. (Errors and 
corrections are discussed in Section 4.2). 

The data are plotted as functions of a/y (which is defined as ui 
henceforth) versus v in Figs.14, 15 and 16. One usually plots the absorp- 
tion per wavelength al = u versus y when analyzing relaxation effects 
associated with a 3; however, accurate velocity data were lacking for each 
frequency at which absorption measurements were made. If it is assumed 
that the velocity dispersion is a small effect compared to the relaxation 
phenomena associated with the sound absorption, the ul versus y plot 
should give essentially the same relative functional behavior as the u 
versus y presentation. (Note: The solid curves are drawn through the 
experimental points; the dashed curves will be discussed in Section 5). 

The temperature dependence of the attenuation along the three axes 


for 6, 10 and 14 MHz is shown in Fig. 17. The values along the a-axis at 
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Table 11 
Sound Absorption (cm!) Along the a Axis 

Temp. Frequency (MHz) 

(°K) 

6 10 14 18 22 26 30 34 

250 0.200 0.484 0.861 

240 0.191 0.422 0.774 0.968 

230 oo. 0.99] 0.702 0.832 0.861 

220 Be 0,356 0.023 0.723 0.821 0.932 

210 0.169 0.343 0.639 0.662 0.750 DON 

200 0,152 02313 0.630 0.576 00035 0.192 

190 Rr 07272 0,563 0.596 0.593 0.721 0.725 0.774 
180 IAS 0,263 0.521 0.536 0.578 0.679 0.625 0.681 
170 ERRO 0,228 — 0.455 0,458 0.458 0.595 0.981 
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Table III 


Sound Absorption (cm!) Along the 
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Table IV 


Frequency (MHz) 


18 


0.806 


0.776 


0.742 


0.675 


0.636 


0.575 


0.939 


0.493 


22 


DR 


02325 


0.699 


0.631 


0.621 


0.619 


02337 


26 


0.994 


0.808 


007593 


0,705 


0.668 


00527 


Sound Absorption (em7!) Along the c Axis 


30 34 38 


0:990 


0.777 0.805 


0.655 0:776 0.930 





| 7 end! 


ZHN)AOINANDIUS 
OS Y20£ 92 22 ai di ol 98 









Zr. 


T »=»0o021! 

qe Mo 06 | 

»40DES< 
SIXD D 

i MoO 
INIZNIS 


AIN3INDIYZW SA ¡7 


of 


> d 


(93S ,.W9) 0601 x ‚TI 


OS 


R 
i v 


ce 


SIXD, q 
3N3ZN38 
AONSNOAYS SA 2 


e 
cy ound 


(ERIN) ADN 
02923 82 et gi el 9 


Lt] 
Et 
Ey 
Ly LL. 
LE 
ES 





=" b è è 9 ` 


H 1002] 
>% 006G] 


AeDle 
A ek é 


ove 





Sa E x © 


TER ees Pag Or ee eT Oe So em t. 


x 


(935- 9) ol 





oj ənd! 
(ZHI) AONSNO ZN J 





OZ! OS sagveossz zz si di 
f A No E 
SE 
AN 
A N 
a Se 
e N 
A N 
v2 a EN 
Ze 70-0 N 





> U x 








¡rl 


x 


(03S WO) 0] 





0.9 


0.87 


97 


0.6 


2.5 


CM 


0.4 


Se 


0.2 


O. I 


2772 


A a vs. TEMPERATURE 
D BENZENE 

| \ O a axis S56. Vie 

IR A b axis OA a 

A S x c axis —- 14 MHz 


yA NS 
290 240 230 220 210 200 190 180 170 


MENCERATURE (°K) 
Figure I7 





Ber 


14 MHz for 200, 210, and 220°K display an anomalous behavior. a remains 
almost constant at these temperatures whereas for all the other measurements 
(a, b, and c axes) a decreases monatomically with decreasing temperature 
at a given frequency in the same temperature range. Whether the anomalous 
effect is real or not is questionable. Since the remainder of the data 
tends to contradict it, it is believed that the values of a along the 
a-axis for 14 MHz at 200 and 210°K may be too high. The results will be 


discussed in detail in Section 5. 
4.2 Errors and Corrections 


In addition to the normal errors due to the experimental apparatus and 
the scatter of the data, there are a number of loss mechanisms peculiar to 
the use of ultrasonic techniques which cause the measured sound absorption 
to be greater than the actual value for the material under consideration. 
Included among these are losses due to diffraction, transducer bonding, 
honparallelism between the sample faces, transducer thickness, reflection, 
and electromechanical conversion (insertion loss). The relative importance 
of each of the effects depends primarily on the magnitude of the intrinsic 
absorption of the sample, the frequency range of the experiment, and the 
Specific acoustic impedances of the transducer material and the sample. 

For the case of solids where sound absorption measurements are made as a 
function of temperature, magnetic field, or some other variably parameter, 
an accurate estimate of all the loss mechanisms is difficult to obtain. 
Consequently, the relative change in a as a function of the parameter of 
interest is obtained and an approximate estimate of the magnitude of each 


loss mechanism is calculated or found experimentally. The errors and 
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corrections pertinent to this experiment will be discussed below. 

The error associated with the experimental equipment depends on a 
knowledge of the sample thickness and how accurately the difference in the 
magnitude of two peaks in each echo pattern can be measured. The sample 
thickness was red ulch a micrometer with within +0.02 cm and was cor- 
rected for thermal contraction by using the values of the linear expansion 
coefficients given by Cox et SID The attenuation could be measured to 
within +0.1 db below 30 MHz and to within +1.0 db at 30 MHz and above. The 
accuracy of the measurements for 30 MHz and above was limited by the low 
signal-to-noise ratio of the patterns obtained at these frequencies. For 
the lowest values of attenuation measured, i.e, at 6 MHz and 170°K, the 
error was calculated to be 4.7% of the measured values. For 30 MHz and above, 
the error estimated is 7.4%. Hence, the systematic error is less than or 
equal to the 7% error associated with the scatter of the data; consequently, 
the latter quantity is the predominant error of the measurement. 

The largest correction--that for diffraction--was calculated by using 
the method of Seki, Granato, and Truell°* as modified for highly attenuating 
Materials by Papadakis.º” Because of the high absorption of benzene, the 


diffraction correction was significant only for the measurements taken at 


6, 10, and 14 MHz. The maximum corrections for these three points were 


29%, 10.3%, and 5.7% of the measured values, respectively (b-axis 1709k): 


At 18 MHz and above the correction was less than 2%. 
Although an estimate of the bond loss was obtained experimentally, an 
idea of the effect of the mylar interface can be obtained by considering 
i : e DOS e 
the standard three-layered-media acoustic transmission problem. fhe bond 


in this casé can be divided into three three-layered-media combinations: 
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(1) transducer-Nonaq-mylar, (2) Nonaq-mylar-glycerine, and (3) mylar- 
glycerine-sample. If a good bond has been made, the Nonaq and glycerine 
interfaces should be thinner than the mylar (less than 1/4-mil). In 
general, it is found that the effect of the intermediate layer depends 
on the magnitude of the propagation vector of the sound, k = 2n/A , and 
m, the thickness of the layer. If the quantity k2 << 1 , the interme- 
diate layer is acoustically transparent. For the mylar layer, the maximum 
value of k& was 1 x 1072 for the frequencies used in this experiment. 
Consequently, if the nonaq and glycerine bonds provide good coupling of 
the sound energy into and out of the mylar seal, the latter should not have 
appreciably interfered with the passage of the sound wave between the 
transducer and the sample. 

To determine the actual effect of the bond, measurements were made on 
a sample which was first air backed (sealed with clear mylar) and then 
backed with a dummy transducer that was bonded onto the reflecting surface 
of the crystal with the same glycerine-mylar-Nonaq bond used for the active 
transducer. The latter configuration yielded higher attenuation measure- 
ments; the excess was attributed to the bond loss. The dummy transducer 
had a fundamental frequency of 5MHz and was one inch in diameter. A trans- 
ducer with a larger diameter than that of the sample was used to eliminate 
the problem of having to center the transducer in the middle of the sample 
so that it would be in line with the active transducer. Alignment would 
be necessary with a transducer which had a smaller diameter than that of 
the sample to ensure that the sound beam is coupled into the dummy trans- 
ducer and not reflected off the benzene surface. 


Measurements were taken between 230-250°K on a nonoriented sample in 
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the two configurations described. The excess absorption due to the dummy 
transducer was approximately 10%. Hence, the bond losses were not a major 


part of the apparent attenuation and the relative behavior of the sound 


_ absorption with temperature should be that characteristic of the sample 


and not that of the bond. 

In general, losses due to nonparallelism of the sample faces are sig- 
nificant at high frequencies and for materials with low intrinsic attenua- 
tion. The frequency at which nonparallelism becomes a predominant loss 
mechanism must be calculated for each case. For highly attenuating materials 
an estimate of the losses due to nonparallelism is given Be 


DOO AD 


s2L 


ol is the attenuation in db/cm due to nonparallelism, R the radius 


of the active area of the transducer, 6 the amount of nonparallelism 


' between the two sample faces expressed in radians, n the echo farthest 


out in the echo train used for measurement, s the sound velocity, and 

L the sample thickness. The value of 6 for this experiment was estimated 
to be 6.67 x 107+ radians. Substitution in Eq. 4.1 indicates that the at- 
tenuation due to nonparallelism at 38 MHz at 180°K is 1.5% of the measured 
value. This represents the maximum error due to nonparal.lelism. 

It was mentioned above that E results were obtained when 
transducers with different fundamental frequencies were used. The same 
effect was observed by ees who made sound absorption measure- 
ments on three different types of materials with five different fundamental- 
frequency transducers. A set of measurements was taken on each material 


with each transducer. For a given type of material, the frequency dependence 
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of the absorption was the same for the five transducers but there was a 
significant difference in the magnitude of the apparent attenuations, the 
thickest transducer (lowest fundamental frequency) yielded the highest 
value of attenuation. For example, the absorption value measured at the 
third harmonic of a 5 MHz transducer was higher than that obtained with a 
15 MHz transducer operated at its fundamental. 

Because the "thickness" effect becomes more predominant with thicker 
transducers, it was a significant factor in this experiment. The 2.0 MHz 
transducer is much thicker than transducers usually employed in making 
sound absorption measurements in solids (normally transducers with canoe 
mental frequencies of 5 MHz and higher are used). An attempt was made to 
obtain an estimate of the thickness effect by applying Sylwestrowicz's 
procedures to this experiment. Measurements were taken with 2.0, 2.5, 
and 5.0 MHz transducers on a benzene sample at 250°K. The details of the 
procedure and the experimental results are given in Appendix G; a possible 
cause for the effect is also discussed in that appendix. 

The experimental results indicated that with a 2.0 MHz transducer, 
the thickness effect accounted for about 30% of the apparent attenuation. 
The experimental values quoted in Section 4.1 were not corrected for this 
effect, however, because the "thickness effect" experiment was carried out 
on one axis at one temperature. The fact that such an effect appears to 
exist is brought out to emphasize the relative nature of sound absorption 
measurements in solids as a function of temperature. One must be careful 
when comparing measurements of different investigators to keep the thickness 
effect in mind, especially when transducers with fundamental frequencies 


below 5 MHz are used. 
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Finally, mention should be made of losses due to reflection and 
electromechanical conversion. For the case of a solid with an air-backed 
reflecting surface, the mismatch of the acoustical impedances between the 
sample material and the air is so great that a reflection coefficient of 
unity can be assumed, er, there are no losses due to reflection. The 
reflecting surface of the sample in this experiment was covered with mylar; 


however, the mylar was so thin that any sound energy absorbed in it would 


| be negligible compared to that of the benzene. Thus, it was assumed that 


the sample was essentially air-backed and the reflection losses were 
negligible. 

Attenuation due to insertion loss can also be neglected in this case. 
This loss can be expressed as the amount of acoustical energy extracted 
from the sample per reflection that is converted into electric energy and 


depends primarily on the electromechanical coupling coefficient of the 


_ transducer employed. For a quartz transducer operating under the conditions 


a 


of this experiment the insertion loss was calculated to be 0.06%. 

Thus, the predominant error in the experimental measurements resulted 
from the scatter of the data with additional uncertainty for the low fre- 
quency measurements (6 and 10 MHz) introduced by diffraction and for the 


high frequency values by nonparallelism (1.5% at 38 MHz) and the low signal- 


- to-noise ratio of the patterns obtained (30 MHz and above). The bond loss 


was found not to predominate over the intrinsic attenuation and the "thick- 
ness effect" was shown to influence only the magnitude of the measurements 
and not the relative frequency behavior with temperature. It is estimated 
that the total error of the measurements is approximately 7% at 26 MHz and 


below and nat greater than 10% for 30 MHz and above. 





ot 
5. Baseussion of Results 
5.1 General 


The sound absorption measurements obtained in this experiment were 
recorded at higher frequencies and lower temperatures then previous mea- 


surements in solid benzene, and, more significantly, were made along known 


| crystallographic directions. A different relaxation behavior was observed 


along each axis which indicates that the sound absorption depends on the 
orientation of the molecules in the unit cell. The lack of relaxation 
effects in the same frequency range in the liquid state further indicates 
that the sound absorption mechanism must be related to the structure of 
the solid, i.e., the ultrasonic beam induces the lattice vibrations to 
interact with each other or causes interactions between the lattice and 
internal oscillations. If a were due to intermolecular vibrations only, 
a reasonable agreement with the theory of Woodruff and Ehrenreich > would 
be expected. However, as pointed out in Section 2.3, Woodruff and Ehren- 
reich's theory predicts the magnitude of a to be less than 1% of the 
observed value and transition times which are three to four orders of mag- 
nitude faster (~ 10-13 sec) than observed in this experiment (» 107º sec). 
Hence, Liebermann's postulate that sound absorption in molecular solids 


arises from the interaction between the inter- and intramolecular vibrational 


‘modes would appear to be appropriate. 


The large values of a observed in this experiment are consistent 
with the magnitude of the sound absorption observed in previous ultrasonic 
attenuation investigations in organic molecules. However, the frequency 


dependence of a above 14 MHz and its temperature dependence below 250°K 
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depart significantly from what has been observed experimentally and from 
that predicted by both the Liebermann and Danielmeyer analyses. In all 
previous investigations, with the exception of the most recent by Wilson 
and Ge, in naphthalene (discussed later), a exhibited an approximate 
frequency-squared dependence and was practically independent of temperature 
within the limits of the experimental measurements. The upper frequency 


and lower temperature limit of these measurements were as follows: 


Benzene” 10.7 MHz 2150K 
Cyclohexane” 8 MHz 185°K 
Carbon ee T359 MHZ LL SS 
p-dichlorobenzene” 65 MHz 220K 
Naphthalene” 95 MHz 2002 


In general, the values of a obtained in this work at 14 MHz and below 
show a partial agreement with the v? dependence of previous investigations; 
however, above 14 MHz, relaxation effects predominate. Therefore, a dis- 
cussion of the data at 14 MHz and below will be given first; then the relax- 
ation effects observed along each axis will be discussed and finally, an 
attempt will be made to relate the behavior of a to the molecular struc- 


ture within the unit cell of benzene. 


5.2 Sound Absorption Behavior Between 6 MHz and 14 MHz 


The temperature dependence of a for each axis at 6, 10, and 14 MHz 
is shown in Fig. 17. In general, the change in a with temperature is the 
greatest for the b-axis and the least for the c-axis but the rate of change 
for each axis increases with frequency. At 6 MHz the value of a for the 


c-axis changes by 4.4% between 215°K (the lower temperature limit of 
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Liebermann's experiment) and 250°K. This compares favorably with Lieber- 
mann's result of a exhibiting only a slight temperature dependence; however, 
the remainder of the data in Fig. 17 contradicts the temperature dependence 
observed in all previous investigations (except for the most recent values 
for E alene”), In the latter case, the values of a also decrease 
with temperature down to the point where relaxation effects are observed 
(170°K) but the rate of decrease does not change as rapidly with increasing 
frequency as it does in the present work. 

It can be seen from Fig. 17 that the relative magnitude of the atten- 
uation is the greatest for the c-axis and the least for the b-axis at 6 
and 10 MHz. The same relative relation holds for the magnitude of the 
linear expansion coefficients B' along each of the axes, i.e., the lar- 
gest value of Bi is that for the c-axis (2.21 x 1074 /0¢)°9 and the smallest 
for the b-axis (1.06 x Ko d a ratio of 2.09. Between 180°K and 250°K 
the ratio of a for the c-axis to that of the b-axis ranges from 1.5 to 
| 2.7 at 6 MHz and 1.1 to 2.3 at 10 MHz. Although DI is not the only 
parameter affecting the sound attenuation, the linear expansion coefficient 
appears to give a qualitative relationship of the relative magnitudes of 
the attenuation along the different axes at 6 and 10 MHz. In general, 
linear expansion in a solid results from the anharmonicity of the potential 
energy which describes the binding forces of the atoms or molecules in the 
substance under consideration. For molecular solids a would also appear 
to depend on the anharmonic terms of the intermolecular potential, conse- 
quently, some correlation between a and DI would be expected, i.e, the 
molecular arrangement within the unit cell does influence the sound absorp- 


tion. At frequencies above 10 MHz the relative relation between a and 
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B' does not hold; however, as will be discussed later, the crystal struc- 
ture influences the sound absorption at these frequencies as well. 

Comparison of the frequency dependence of the absorption values obtained 
at 14 MHz and below with that obtained from Liebermann's values for benzene 
can be. made by Feen the slopes of the lines that result from a loga- 
rithmic plot of a versus v. The slope should be close to two if a 
exhibits approximate frequency-squared dependence. At 250°K the slopes 
obtained from the data of this experiment between 6 and 14 MHz are 1.73, 
1.85, and 1.42 for the a, b, and c axes, respectively. At 170°K the slopes 
obtained from the data in the same frequency range are 1.60 and 1.83, respec- 
tively, and for the c-axis at 180°K, 1.08. Hence, the frequency dependence 
along the b-axis exhibits essentially no change between 170°K and 250°K 
whereas it changes by 7.5% and 24% from the high temperature values for the 
a and c axes, respectively. From Liebermann's measurements on nonoriented 
benzene crystals at 6.4 and 10.7 MHz at 273°K, the slope is 1.90. If it is 
- assumed that the slopes change linearly with temperature, the extrapolated 
Miles for the a, b, and c axes at 273°K become 1.76, 1.85, and 1.52, respec- 
tively. Hence, of all the measurements, the frequency dependence of a 
for the b-axis agrees most favorably with that observed by Liebermann and 
the temperature dependence compares least favorably with Liebermann's re- 
sults, i.e., there is a 62% change in a between 170°K and 250°K (6 MHz). 
| The agreement of the frequency and temperature dependences for the c-axis 
values are reversed, i.e, agreement of the frequency dependence is the worst 
and that of the temperature dependence is the best. 

The fact that the slopes of the lines on the logavs, log v plots for 


the a. and: ec axes decrease with decreasing temperature suggests that 
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presence of relaxation effects; this is borne out by the data obtained 
above 14 MHz (see Figs.14-16). Wilson and Yun observed the same behavior 
for naphthalene; i.e., the slope changed from 1.74 at room temperature to 
approximately one in the temperature range where relaxation effects were 


observed (170°K). 
5.3 Pertinent Facts Concerning Relaxation Theory 


Prior to discussing the experimental results for the higher frequencies, 
a brief description of the standard relaxation curve on a u (absorption 
per wavelength) versus v plot will now be given. The standard formula 


for the attenuation per wavelength (Eq. 2.37) can be rewritten in the form 


TY Um 
u = CD A (5.1) 
(1-r) Ltrs 


/2 


where T= Tas (=) . Thus, u can be separated into two terms, one 
which involves only the relaxation strength r and the other which involves 
wT The latter term controls the shape of the curve on the y versus v 
plot. If the quantity wt /(tw*t?) is plotted as a function of wt» a 
broad bell-shaped curve is generated which has a maximum value at 

ut = 1. In the region where SS the curve is a straight line with 

a slope of one and where ur >> 1 the curve is also a straight line but 
with a slope of minus one. The curve thus formed is the "standard relaxa- 
tion curve" for a single relaxation process. The position at which wt = 1 
on the u versus v plot defines the maximum absorption per wavelength Um 
and the relaxation frequency Vo? The change of Um and vo with temper- 


agure can be used to obtain information about the reaction mechanisms which 
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take place in the relaxation region. If velocity dispersion is small, y 
can be replaced by u' , i.e. (a/v) , and the u! versus v plot should 
exhibit the same characteristics as the u vs. v presentation. 

A single relaxation process will occur if there is a maximum value of 
sound absorption per wavelength that results from the exchange of energy 
between one lattice vibrational mode and one optical vibrational mode. As 
discussed in Section 2.2, in many cases in liquids and gases the entire 
optical spectrum of a substance will relax as a group; such a relaxation 
will also generate a standard relaxation curve on the u vs. v plot. If 
there is more than one relaxation process occurring at frequencies which 
are relatively close together, i.e, which differ by less than a factor 
of about 10, the standard relaxation curve will become distorted. One can 
‚, usually obtain a reasonable reproduction of the distorted curve by adding 
two or more standard relaxation curves with different values of Um and 
E However, multiple relaxation processes are quite complicated and the 
addition procedure only gives an indication that such a process exists; it 
does not yield any information as to how the interaction is taking place; 
i.e, whether the processes are interdependent or independent of each other. 
Hence when the additive procedure is used, it must be kept in mind that the 
values of Mn and E of the individual components of the sum describe 


only the simplest possible interaction that is likely to occur. 
2.4 Relaxation Effects in Benzene 


Figure 14 shows the ui vs. v curves for the a-axis. The standard 
relaxation curve has been drawn through the experimental points for 170, 


190, 230, and 240°K. (The trend of the curves at the lower temperatures 
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was used to draw the curves for the upper temperatures.) The curves indi- 
cate that the relaxation frequency vm changes from 17 MHz to 12 MHz and 
mm decreases by 58.5% from its value at 240ºK. The shift of vm toa 
lower value at the lower temperature is characteristic of vibrational relax- 
ations, i.e., normally Te transition times will increase with decreasing 
temperature. 

All the points do not fall on the curves as drawn, however. The data 


from 18 MHz to 34 MHz fit very well as does the 6 MHz point, but the 10 MHz 


and 14 MHz points are below and above the curves, respectively. The 


discrepancy at 14 MHz becomes more pronounced as the temperature is decreased. 


A better fit to the points would be given by a curve which peaks more 


sharply than the standard relaxation curve; the dashed curve drawn through 


A 


| the points for 190°K represents such a behavior. 


The experimental points for 190°K have been replotted in Fig. 18. The 
curve drawn was obtained by adding the contributions of the three standard 
relaxation curves shown and appears to fit the experimental data somewhat 
better than the standard curve. This suggests the possibility of multiple 
relaxations occurring in the same frequency region. For the example shown, 
the relaxation frequencies would be 9, 14, and 24 MHz and the relaxation 


strength of the center curve is about 80% higher than the other two. As 


indicated earlier, this is only representative of the type of multiple 


relaxations which might occur. 

The data for the b-axis are shown in Fig. 15. Every experimental point 
except one (10 MHz @ 170°K) falls on the standard relaxation curves drawn. 
Once again the trend of the low temperature data was used to draw the curves 


for the high temperature data. The points at 170°K fall on the back side 
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of the standard curve which gives a good indication that a relaxation.is 


| present. The curves indicate a relaxation frequency of 26 MHz at 240°K 


and 20 MHz at 170°K with a change of 53.5% in un over the same temperature 
range. The shape of the curves indicate the presence of a single relaxation 
along this axis, i.e., either one lattice mode is exchanging energy with 

one optical mode or the entire optical spectrum is relaxing as a group with 
a single relaxation time. 

The c-axis curves shown in Fig. 16 exhibit the most distortion of all 
the curves obtained, especially at the lower temperatures. The experimental 
points for 240°K roughly fit a standard curve centered at 14 MHz; however, 
as the temperature is lowered another relaxation appears to influence the 
absorption. At 180°K the points from 26 to 38 MHz appear to be on the up- 
ward slope of a second relaxation. In Fig. 16, two standard relaxation 
curves have been drawn in below the experimental curves. By adding the 
components of these two curves, the resultant curve (see Fig. 19) qualita- 
tively follows the experimental points at 180°K. This suggests that there 
is a low frequency relaxation at 4.5 MHz with u! = 3.5 x 1078 em! sec 
and a high frequency relaxation with a relaxation frequency of 120 MHz and 
un = 3.0 x 10°® cm”! sec. Once again, the behavior of E follows the 
pattern which is characteristic for vibrational relaxations. At 240°K the 


high frequency relaxation peak is at a value of ve and um such that the 


values of ui associated with that curve do not influence the values of 


ui of the relaxation curve centered at 14 MHz (within the range of the 
experimental measurements). As the temperature is lowered both relaxation 
peaks shift to a lower frequency. The influence of the upper peak starts 


to appear at 210°K and becomes stronger as the temperature is lowered further. 
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The exact behavior of the upper relaxation cannot be determined because of 
the limited experimental information. However, the low frequency relaxation 
appears to be centered at 14 MHz at 240°K and at 4.5 MHz at 180°K; ui 

m 


changes by 32.8% over the same temperature range. 
5.5 Comparison With Theory 


Neither Liebermann's nor Danielmeyer's theory predicted the existance 


| of relaxation effects in the frequency and temperature ranges of this ex- 


periment. In Section 2 it was pointed out that both theories contained a 


number of fundamental weaknesses and didnot adequately describe the sound 


absorption process in molecular crystals; however, since they are the only 


theoretical treatments available, a brief comparison will be made between 


' each theory and the experimental results. In so doing, modifications will 


' be made in the theories to take into account the anisotropic properties of 


the solid wherever possible and factors in the equations for a and fT 
that may be treated as experimentally adjustable parameters will be pointed 
out. 

A comparison of the experimental relaxation times, i.e., where wT = » 
for each axis and those calculated from Liebermann's and Danielmeyer's 
theories is shown in Table V. The comparison is made for temperatures of 


180°K and 240°K, the temperature range which is common to the data for all 


three axes and for which an experimental value of Wi could be obtained for 


each axis. (Fewer points were obtained at 250°K than at 240°K and below; 
consequently, a value of va could not be determined for each axis at 
250°K). Equations 2.5 and 2.26 were used to calculate 1 for Liebermann's 


and Danielmeyer's theory, respectively Lu, was considered to be constant 
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Table V 


Comparison of Experimental and Theoretical Relaxation Times 


Temp. ES vi tia T 
(°K) (sec x 10°) 
Experimental 240 9239. 26,15. 21124 


180 13.4 8.00 392) 


Danielmeyer's 240 02267 2.47 10 103 
Theory 180 0.45 4.74 ¿0073 
Liebermann's 240 0.94 
Theory 180 3.70 


for the present case). Since Danielmeyer's theory is derived from a solid 


state viewpoint, it is possible to modify it slightly to account for aniso- 


tropy. To calculate t from Eq. 2.26 one must first calculate the value 
of the coupling parameter B contained in the factor K (see Eq. 2.23). 

B is related to the Griineisen constant by Eqs. 2.20 and 2.21. Hence, the 
Griineisen constant can be calculated as a function of temperature for each 
principal direction from Eq. 2.21. (In doing this 8 must be replaced by 
38' for each direction and e by the experimentally determined values 
of the bulk modulus (Heseltine et ADE Once this is done, Eq. 2.20 can 


be used to calculate the value of B as a function of temperature using 


- the value of Be which is appropriate for a longitudinal sound wave propa- 


gated along the axis under consideration (elastic constants also given by 


Meseltine et alt"). Since the coupling parameter A contained in Lieber- 
mann's equation for t (Eq. 2.5) is derived from the Lennard-Jones potential 


a quantitative relation between A and the elastic and thermodynamic 
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properties of the crystal cannot be made; consequently, the value of A 
originally calculated for benzene by Liebermann was used. The relaxation 
time predicted by the modified Danielmeyer theory for the b-axis agrees 
fairly well with the experimental values for the same axis. Agreement for 
the a and c axes is much poorer. The value of t predicted by Lieber- 
mann's theory agrees to within one order of magnitude with all the experi- 
mental values. Except for I both theories predict an increase in the 
relaxation time with decreasing temperature as is indicated by experiment. 
A indicates the opposite behavior and results from the fact that the 
values of K for the c-axis decrease much more slowly with temperature than 
for the other two axes. This allows the exponential factor Se (see 


Eq. 2.29) to have a greater influence on the change of 1 with decreasing 


' temperature. 


Calculated for the a-axis is divided by two, the predicted values of 1 


Better agreement between the value of ER and To obtained from 
experiment and those predicted by Danielmeyer's theory can be obtained by 
treating B as an adjustable parameter. In Eq. 2.20 only Eh and E 
were used to calculate the values of B for the a and c axes, respec- 
tively; however, as pointed out in Section 2.2, B is related to the second 
and third order elastic constants. For an anisotropic solid the value of 
B may be quite different from that given by Eq. 2.20. If the value of B 
a 
fall into the same order of magnitude as the experimental values. For the 
c-axis, the value of B for 240°K and 180°K must be divided by about 3 and 
3.5, respectively so that the predicted values of ie will be of the same 


order of magnitude as the observed values and also increase with decreasing 


temperature 
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The above discussion assumed Vo to be constant; however, as mentioned 
in Section 2, when Yo is treated as a function of temperature, Liebermann's 
theory predicts the attentuation to decrease with increasing temperature and 
Danielmeyer's theory predicts an increase in the relaxation time with in- 
creasing temperature (increase of a with increasing temperature). Although 
Danielmeyer's results is plausible, it does not agree with the relaxation 


behavior observed. For a decrease in the relaxation time with increasing 


temperature would shift the relaxation peak to a lower frequency as the tem- 


perature is increased, which is just opposite the normal behavior in vibra- 


tional relaxations and from that experimentally observed. 
Another parameter which can be compared with the theory is the change 


in He for each direction as a function of temperature. To do this the 


relaxation strength in Eq. 5.1 must be modified to take the sound velocity 


into account. In general H occurs when ot = i= and For this Case Eq. 5,1 


becomes 


Tr 
Ep m... 2), 
E EE 


From Danielmeyer's theory the relaxation strength is given by Eq. 2.38 and 
for benzene yields values of r << 1 (% 1073); therefore, the r in the 
denominator of Eq. 5.2 can be neglected. Equation-5.2 can be rewritten in 


terms of D as 


t z £ AA 5.3) 
Ve 2 (=) S E ( 


Here s is the longitudinal sound velocity along the principal axis for 
which the absorption measurements are being made. 


Equation 2.38 shows that r is temperature dependent, i.e., the specific 


heats and the factor z_ , which is related to the elastic constants 
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(Eq. 2.39), change with temperature. Hence, H, can be calculated for 


each direction at different temperatures which enables a comparison of the 


change in the relative values of H. with temperature. 

Since E. is based on an isotropic solid, a slight modification is 
necessary to use it in the present case. The principal axes of the ortho- 
rhombic structure are two-fold symmetric; hence, there are two shear veloci- 
ties for each principal direction. To calculate E for a given direction 
for benzene, the average value of the two shear velocities was used. Equa- 


tion 2.39 becomes 


(5.3) 


Here ec and cm are the elastic constants associated with the shear 
1 2 


velocities which propagate along one of the principal axes; > is the 


corresponding elastic constant for longitudinal waves. Hence, to calculate 


D, Eq. 2.37 is modified to become 


i 
Ce ae 
[Mees ert = — O GE A (5.4) 
_rtaty _ _ytyot 
Ee. IG? Gell GG, 


RE e 
By using the experimental values of the sound velocities and specific 


a Um can be calculated from Eq. 5.4 for any temperature between 


heats 
HECK and 250°K. 
A comparison of the change in H, between 180°K and 240°K that was 


experimentally observed and that calculated from Eq. 5.4 is given in Table VI. 


It can be seen that the predicted change in un is from 1.55 to 2.12 times 


greater than that obtained from the observed values. A better agreement 


could be obtained in two ways. First, the factor be can be treated as an 
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Table VI 


Comparison of the Theoretical and Experimental 


Change in un With Temperature 


Temp 
(°K) 
u (exp) 240 
(cm-!sec x 107) 180 
q 
y, (th) 240 
(cm~!see x 107) 180 
q 
Ay (exp) 
q 
Ay, (th) 
Ay, (exp) 
Ratio 


1 
Ay, (th) 


4.40 


122 


uu. us 


71252 


To! 


1.55 


Se 


EG) 


32.0% 


68 % 


2.12 
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experimentally adjustable parameter since it is probably related to the 
elastic constants of the material in a more complicated way aa is described 
by Eq. 5.3. It would be possible to adjust ce For each axis and for each 
temperature to obtain the change in De that is experimentally observed. 
However, this procedure does not yield any information concerning the details 
of the physical processes involved. 

A second way to obtain better agreement is to consider the possibility 
of partial relaxation effects; i.e., the energy that is coupled into the 
internal modes does not relax with all the modes at one time as assumed by 
Danielmeyer. This type of behavior has been experimentally observed for 
vibrational relaxations in liquid sulfur dioxide (Bass and Lamb?) and 
liquid methylene chloride (Andreae°”), 

For liquids and gases the relaxation strength and the internal specific 
heats are related by an expression similar to Eq. pees However, rather 
than use the total internal specific heat, the contribution of the internal 
specific heat due to each mode is determined from the Einstein relationship 
mq. 2.3) and D, is calculated for the mode under consideration from the 
result. Absolute sound absorption measurements are usually made in liquids; 
consequently, it can be determined if the total internal specific heat 
relaxes at the experimentally observed relaxation time by comparing the 
numerical value of Um obtained from experiment to that calculated from 
theory. For liquid sulfur dioxide and methylene chloride it was found that 
u (exp) was less than u (th) when the latter was calculated by use of the 
total internal specific heat. However, if only the contributions to si 
due to certain internal modes were used to calculate u (th), excellent 


agreement between theory and experiment was obtained. For sulfur dioxide 
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there are four internal modes (one of which is degenerate) and it was found 
that only two contribute to the observed relaxation. For methylene chloride 


only eight of the nine internal modes participate in the observed relaxation. 
5.6 Relationship of a to the Molecular Vibrations 


A procedure similar to that employed for liquids could be used to cal- 
culate ut (th) in the present case by replacing C; in Eq. 5.4 by the 
C; (calculated from Eq. 2.3) for the internal modes of crystalline ben- 
zene. A comparison between y, (th) calculated in this manner and u' (exp) 
could be made on a relative basis for two different temperatures as was done 
in Table VI. Before going into the details of this calculation, a brief 
description of the benzene molecule will be given. 

Cox et a determined that the benzene molecule is accurately planar 
with one carbon atom at each vertex of a regular hexagon. One hydrogen atom 
is bonded to each carbon atom and the C-H bonds extend radially outward from 
the center of the hexagon along lines joining the center and the vertices 
of the hexagon. A schematic diagram of the molecule is shown below. 


H, 





H4 


Fig. 20. Benzene Molecule 
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The numbering system is that used by Cox to describe the position of the 
atoms when the molecule is in the unit cell E crystalline benzene (dis- 
cussed later). The internal oscillations of the benzene molecule can be 
divided into in-plane and out-of-plane vibrations, the latter usually cor- 
respond to the lower EE EE es in the optical spectrum; hence, are the 
ones with weaker bonds. Examples of the out-of-plane vibrations would be 
the C-H out-of-plane bending mode and the C-C-C-C skeletal out-of-plane mode. 
The weakest in-plane vibration is the C-C-C bending mode which tends to 
change the obtuse angles of the type C,-C,-C,. The strongest bond of all 
the internal vibrations, i.e., that with the highest frequency, is the C-H 
in-plane stretching vibration. There are a total of thirty normal modes 
for the internal oscillations of the benzene molecule (for a complete de- 
scription of each mode see Wilson et ai, Dn. In the gaseous and liquid 
states ten of the modes are e however, in the solid state the 
degeneracies are lifted and there are thirty nondegenerate internal nales 
Many spectroscopic experiments have been performed on crystalline ben- 
zene; consequently, the internal vibrational frequencies are well known. 
By using the most recent Raman ‘+ and EE data for benzene SC. was 
calculated for the internal modes in the lower part of the optical spectrum 
and for the highest internal mode. Then GH was calculated using the 
Definition in Eq. 5.4 except that ça was replaced by E 
ôC.. From this a value of y, (th) can be calculated for each internal mode. 
The results are shown in Table VII. The modes are listed in order of in- 
creasing wavenumber, the value in parentheses corresponds to the standard 


S S 70 
numbering system for the normal modes first introduced by Wilson. Also, 


the type of ‘vibration is indicated; i.e, C-H bending, etc. 
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It can be seen in general that the change in D, between 180°K and 
240°K becomes greater as the internal vibrational frequency increases. For 
the a and b axes the changes in Um for the lowest optical mode are 
within a few percent of the experimentally observed values. The correlation 
for the c-axis is not as good; however, it is much better than that obtained 
in Table VI. Hence, Table VII gives an indication that the lower frequency 
internal modes are the ones which play the principal role in the relaxations 
which are observed. Also, it should be noted that the values of the y (th) 
calculated with the 6C are of the same order of magnitude as the ut Cexp) 
whereas the values of u (th) calculated with the total internal specific 
heat are one order of magnitude greater than ut (exp) (compare Tables VI 
and VII). Although it was indicated in Section .2 that the experimental 
values of a could be as much as 40% higher than the absolute values, i.e., 
due to the thickness effect and experimental error, the fact that the mea- 
sured values are significantly less than those predicted by using the total 
internal specific heat gives further support to the idea of partial relax- 


ation processes. 
5.7 Relationship of a to Crystal Structure 


An intuitive argument as to why one might expect the low frequency 
internal modes to participate in the relaxation processes rather than the 
other modes can be obtained by trying to visualize the effect of the sound 
waves on the molecules in the benzene unit cell. The arrangement of the 
benzene molecules in the unit cell as seen looking down the b and c axes 
is shown in Fig. 2l. The crystal structure of benzene is simple ortho- 


rhombic with a space group of Pbca which means there is a molecule 








Fig. 2la. Diagram of the crystal structure of benzene 
viewed down the b-axis. 





Fig. 21b. Diagram of the crystal structure of benzene 
viewed down the c-axis, 


Fig. 21. Crystal structure of benzene (from Cox, Cruikshank, 
and Smith [60]). 
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positioned at each corner and in the center of each face of the unit cell. 
(The structure can be thought of in simpler terms as face centered ortho- 
rhombic). Hence, there are four molecules per unit cell; each molecule has 
twelve nearest neighbors. The four molecules have been designated I, II, 
III, and IV in Fig. 2la. The molecular planes are almost parallel to the 
b-axis and molecules of type I and II are inclined about 45º to the a and 

c axes and are at right angles to each other. (The same relative orienta- 
tion exists for molecules of type III and IV.) The He hydrogen atom (see 
Fig. 21b) of molecule I is located between the H, and H. hydrogen atoms of 
molecule II; consequently, the C-H bonds of the two molecules intermesh in 

a bevel gear arrangement. The bevel gear arrangement forms a line along 

the a-axis in the shape of a corrugated sheet. Molecules III and IV form 

a Similar sheet which appears to be stacked on top of the sheet formed by 

the type I and II molecules. Thus, there are a series of corrugated sheets 
perpendicular to the c-axis. 

Lattice vibrations in a molecular solid can be divided into transla- 
tional and rotational modes. In benzene the rotational modes are the lattice 
vibrations with the highest KEE thus, are the ones in the upper 
Bert Of the acoustic spectrum. Cox et a have determined that at thermal 
equilibrium the molecules rotate about mutually perpendicular axes in the 
molecular plane and about the sixfold axis perpendicular to the molecular 
plane. The latter rotation has the greatest amplitude of all the rotational 
modes and causes the type I and II molecules to oscillate back and forth in 
Opposite directions; i.e., as two gears would if one of the gears were set 


into simple harmonic motion. This action causes all the molecules in the 


same sheet to oscillate in the same manner. 





SES 


If a longitudinal sound wave is propagated along the c-axis (see 
Fig. 2la), the ensuing compressions and rarefications will tend to induce 
the molecular planes to oscillate about axes which are parallel to the 
b-axis. Thus, the corrugated sheets will be flattened out in the regions 
of compression and will tend to their equilibrium positions in the regions 
of rarefaction. Such a motion would tend to distort the benzene ring and 
also cause the nearest neighbor H-H distances to be displaced from equi- 
librium. One would expect that the internal modes most affected would be 
the skeletal out-of-plane vibrations, the C-C-C in-plane vibration, and 
the C-H out-of-plane bending vibrations. From Table VII it can be seen 
that these three modes are Ene lowest in he optical spectrum. Thus, there 
is a possibility that relaxations are taking place between an excited state 
of the rotational lattice modes and excited states of the lower modes in 
the optical spectrum. 

A similar situation would exist when a longitudinal wave is propagated 
along the a-axis. The molecular planes would be induced to oscillate about 
axes parallel to the b-axis and in regions of compression would be squeezed 


together in an "accordion-like" fashion. Since the interatomic bonds that 


are formed by atoms that are closest together in nearest neighbor molecules, 


D 


4 and H atoms, are much 


e. He atom of molecule I om. and the H 
stronger than those bonds formed by atoms which are farther apart, this 
type of motion would tend to distort the molecular plane (skeletal out- of- 
plane vibration) and induce the C-H bonds to be bent out of the molecular 
plane. Once again, there is the possibility that lattice rotational mode 


is exchanging energy with the lower modes in the optical spectrum. 


A longitudinal sound wave propagated along the b-axis would induce 
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the molecule to oscillate about its sixfold axis. This would cause the 
"gears" in each sheet to oscillate in phase, i.e., the type I and II mole- 
cules would oscillate in opposite directions. Such a motion would tend to 
change the H-H distances between nearest neighbor molecules but the rota- 
tional motion should not cause the molecular ring to be distorted as in the 
previous two cases. However, the molecular ring could be distorted by the 
translational motion of the molecules which would cause the C-C-C in-plane 
vibration to become excited. Thus, the sound wave along the b-axis appears 
to affect the molecules in a much different manner than when propagated 
along the a and c axes. The observed relaxation for the b-axis is also 
different from those of the other two axes. Figure 15 gives the indication 
eat a single relaxation is taking place. This could be a result of the 
interaction of a lattice mode with either a C-H bending mode or the C-C-C 
in-plane mode. It is more likely that the entire optical spectrum is re- 
laxing as a group; however, the comparison in Table VI tends to contradict 
this possibility. If one believes that Eq. 5.4 gives a reasonable value 
for ualth), then a relaxation between one rotational lattice mode and one 
low frequency optical mode is indicated. 

The above discussion, which attempts to relate the action of the sound 
wave on the molecular structure is of course oversimplified. The actual 
interaction in the crystal is certainly much more complicated; however, if 
molecular resonance is the primary sound absorption mechanism in organic 
crystals, one can see that the molecular arrangement within the unit, cell 
should be an important factor in the sound absorption behavior. It would 
be interesting to make sound absorption measurements on an oriented organic 


crystal with cubic structure. In such a case, the elastic constants along 
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all the (100) axes would supposedly be the same and any difference in the 
sound absorption behavior along each principal axis would be a result of 
the molecular arrangenent. 

The agreement between theory and experiment is not particularly good. 
Table V shows an order of magnitude agreement between some of the theo- 
retically predicted and experimentally observed relaxation times; however, 
these values were calculated under the assumption that v, is independent 


D 
of temperature. When the temperature dependence of vn is taken into 
account, Liebermann's theory predicts relaxation times which lead to values 
of a(T) which contraduct experimental results (see page 11). Danielmeyer's 
theory predicts the value of T to decrease with decreasing temperature 
ich is also in opposition to the experimental result. However, both 
relaxation time theories were based on the high-temperature Debye approxi- 
mation, 1.€., 9, << T. The lower temperature limit in this experiment 
(170°K) was close to 8, for benzene (150°K); consequently, the theoretical 


D 
formulas for 1 should not be expected to be valid in that temperature 
region. 

The comparison shown in Table VI suggests that Danielmeyer's assumption 
that the entire optical spectrum relaxes with a single relaxation time is 
questionable. One would expect better agreement between the theoretically 
predicted and experimentally observed change in ua Closer agreement was 
obtained by considering the possibility of partial realxation processes 
(Table VII). The analysis presented in Tables VI and VII depends on the 
validity of Danielmeyer's expression for r (Eq. 2.38). Since Eq. 2.38 


is based on well-known concepts of solid state and relaxation theory, at 


least an order-of-magnitude agreement between r'(exp) and r'(th) should 





-111- 


be expected (using both SR and 6C.). To apply the theory to a specific 
crystal structure, however, a more realistic expression for the factor 
E (Eq. 2.39) must be obtained. By considering an anisotropic solid, such 


an expression could be derived using Danielmeyer's approach. 
5.8 Suggestions for Future Research 


Further theoretical and experimental work concerning sound absorption 
in molecular crystals is desired. The existing theories contain too many 
generalizations to enable accurate predictions of the sound absorption 
behavior. One of the basic problems in the derivation of transition proba- 
bilities was the lack of an accurate expression for the intermolecular po- 
tential. The potential obtained by Williams“ (Eq. 2.8) should be useful 
in this regard. As a first approximation one could take the third derivative 
of Eq. 2.8 (restricting the summation to include only interatomic distances 
of the nearest neighbor molecules (~ 4A°)) and use first-order time-depen- 
dent perturbation theory to calculate tT. Further refinement could be made 


by separating each interatomic distance, i.e., Too? and r 


Pa c CH? into 


components along each axis, then taking the third derivative of each re- 
sulting expression. To calculate t , however, the Debye approximation 
for the density of states would still probably have to be used. The resul- 
tant relaxation time could be related to the complex sound velocity in an 
anisotropic crystal by using an approach similar to that used by Danielmeyer 
(see page 23). 

Further experimentation would be desirable with organic substances 
which have simpler molecular structure and are cubic in the solid state. 


Carbon tetrachloride (CCl, ) fulfills these criteria and large single crystals 
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of ccl, have been BN However, cel, freezes at -23°C and undergoes a 
phase transition at -48%C, The low freezing point would present experi- 
mental difficulties during orientation and the phase transition would limit 
the temperature range of the experiment to about 25°. 

Although a more complex molecule than benzene, cyclohexane forms as 
fcc and its freezing point is 6.2°C. Hence, the crystal growing apparatus 
“and the sample preparation procedure developed for this experiment could 
be used. Orientation would have to be accomplished using X-ray techniques, 
however, since cubic crystals cannot be oriented with a polarizing micro- 
scope (cubic crystals appear to be isotropic when examined with polarized 
light). 

Further absorption and velocity measurements on benzene would be 
desirable. Velocity measurements made as a function of frequency would 
indicate the importance of dispersion in the region where relaxations were 
observed in this work. Absorption measurements at higher frequencies would 
be possible at lower temperatures and may show the existence of relaxation 
phenomena in the 100 MHz region as was qualitatively indicated in this work 
for the c-axis. Shear wave absorption measurements should give further 
indication as to the importance of the molecular arrangement in the unit 
cell. The corrugated-sheet arrangement described earlier indicates that 
molecules of adjacent sheets are much more loosely bound than those within 
the same sheet. This suggests the possibility that the sheets may be able 
to slide with respect to each other; i.e., that the a-b crystallographic 
plane is a slip plane. If dislocation damping is a significant sound ab- 
sorption mechanism in benzene, transverse waves that are polarized in the 


Slip plane should yield higher absorption values than transverse waves 
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polarized perpendicular to it. All the experimental apparatus developed 
for this experiment could be used to make shear wave measurements in ben- 
zene except that an improved transducer bonding technique would have to be 
developed. (Preliminary work with thick shear wave transducers indicated 
that ae bonding technique described in Section 3.6 would not support shear 
waves. ) 

Hence, there are at least four experiments, which would yield useful 
information about the properties of molecular crystals, that could be per- 
formed by using all or part of the apparatus and procedures developed for 
this work: É | 
(1) Sound velocity measurements on oriented benzene crystals in the same 

temperature and frequency range of this experiment. 

(2) Longitudinal-wave sound absorption measurements on oriented benzene 
crystals below 170°K. (Requires development of a new bonding technique.) 

(3) Shear wave sound absorption measurement on oriented benzene crystals 
below 170°K. (Also, requires an improved bonding technique. ) 


(4) Longitudinal wave absorption measurements on oriented cyclohexane 


crystals. (Requires development of an X-ray orientation technique. ) 
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Appendix A 


Procedures Used With the Vacuum Sublimation Apparatus 


To start the sublimation process (see Fig. 1 in Sec. 3.2) all the 


en lon stopcocks are closed except T 


To» and T The storage tube is 


Sp 8° 


put into position by connecting the joint d The system is then evacuated 


to a pressure of 107? mm Hg (lu), which seals Ja: Approximately 200 ml. 


of liquid benzene are put into the supply tube and this tube is attached 


at J The benzene is then slowly frozen by means of liquid NM: Che 


a 9? 


liquid N, is then removed and the stopcock T 


2 is opened. The initial 


5 


surge of air and benzene vapor is trapped by the cold trap immediately in 
front of the rough vacuum pump. When the pressure becomes less than 2 mm Hg 


(20004), the sublimate trap is cooled by liquid No. The level of the liquid 


N, must be maintained slightly below the lower end of the inner tube. If 


the level is higher, the inner tube will very rapidly become clogged and 
the maximum capacity of the trap will not be utilized. When about half of 


the sublimand in the supply tube is gone, stopcocks To and Tg are closed. 


The remaining sublimand is kept frozen with liquid N Liquid N, is then 


93 
removed from the sublimate trap and transferred to the storage tube. The 
resulting pressure difference between the storage tube and the sublimate 


trap causes the sublimate to be "vacuum transferred" into the storage tube. 


When the transfer is complete, stopcock T, is closed and T, and T, are 


7 8 


opened. The liquid N is removed from the supply tube and transferred 


2 
to the sublimate trap; then the remainder of the sublimand is sublimated. 
During this step, the benzene in the storage tube is kept frozen with liquid 
No. When the sublimand is completely gone, the sublimate is again "vacuum 


transferred" to the storage tube. 
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It takes approximately four hours to process 100 ml. of sublimand 
and about the same amount of time to transfer to the storage tube. Since 
the capacity of the crystal-growing jar is approximately 550 ml., it is 
necessary to degas about 600 ml. of liquid benzene. In order to do this, 
three storage tubes were constructed and a different storage tube was at- 
tached to Ja each time the supply tube was refilled. Thus, the above 
procedure is repeated three times. The supply tube, sublimate trap, and 
the storage tube (in position for transfer) are shown in Fig. A-l. A special 
clamp had to be constructed to hold the supply tube in position at Jo 
The storage tube, supply tube, and clamp are shown in Fig. A-2. 

Once the benzene is degassed, it has to be put into the crystal-growing 
jar without allowing air to get into the system (see Fig. 1). Before putting 
the benzene into the crystal grower, the entire system is evacuated, by use 
ef both the diffusion pump and the roughing pump. Stopcocks In and Te 
are opened and the remaining ones left closed. [It was found that the 
lowest pressure attainable with the teflon stopcocks in the system (without 
vacuum grease) was 5 x 107% mm Hg.] After the pressure has reached the 
lowest value, In is closed and the diffusion pump is turned off. A stor- 
age tube is then put into the filling position as shown in Figs. 1 and A-3. 

T, is opened to allow the air between T, and T, to be evacuated and 


3 1 2 


Jy to be sealed. T, is then closed, and T, 


Start the benzene to flow. T, is opened to allow the benzene to be pulled 


is opened just enough to 


into the crystal-growing jar by gravity and pressure differential. After 


the flow is started, T, is fully opened. 


When the transfer is complete (it takes about 1 minute), T and then 


T, are closed and the empty storage tube is removed. Any residual liquid 





-116- 


benzene left between Jy and T, is swabbed with cotton Q-tips; then 
another storage tube is put into position and the transfer process is 
repeated. During transfer with the third storage tube, the flow is con- 


trolled with Ty) and T, to obtain the proper liquid level in the crystal. 


grower as described in Appendix B. 





Supply Tube Storage Tube Sublimate Trap 


Fig. A-1. Vacuum sublimation apparatus 
(Also see Fig. 1). 





“ll 7 = 





Pig. A-2, Storage Tube, Clamp, and Supply Tube 





Fig. A-3. Storage Tube In Filling Position 
(Also see Fig. 1) 
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Appendix B 


Crystal Growth Procedure and Associated Problems 


To start the growth procedure, the jar is filled to a level which covers 
about half the teflon collar. (The crystal growing apparatus and filling 
procedures are described in Section 3.3 and Appendix A, respectively). The 
refrigerator and the two SCR controllers are turned on immediately after 
the Jar is filled. The controllers are set to establish an initial equi- 
librium air temperature of approximately 4°C and a melt temperature of 3.5°C 
at the copper probe tip after the refrigerator has cooled down to its lowest 
temperature (-5°C). (This corresponds to settings of 7.95 and 6.50 on the 
air and melt SCR controllers, respectively.) After the refrigerator has 
been in operation for one hour, the probe and lid temperature Variacs are 
turned on. (This is not done earlier,as a safety precaution, since benzene 
is extremely inflammable and has a low flash point of - 11.1°C.) The Variacs 
are set to establish an initial equilibrium temperature of 1°C on the probe 
and a lid temperature of 8.5°C (settings of 9.8 volts and 20.0 volts respec- 
tively). It takes approximately three hours for all the temperatures to 
come to equilibrium. Thereafter the seed growth is commenced. The melt 
temperature is then decreased, using the motor drive on the SCR controller, 
and a seed should form within twenty-four hours. However, the first seed 
usually is polycrystalline and consequently it has to be partially melted 
by means of the melt heater control. The melt temperature is raised to a 
point where the seed is just barely visible; then the melt temperature is 
again reduced. Usually this procedure will yield one or two single crystal 
on the probe tip but three or four eee may be necessary to obtain one 


that is satisfactory. After the seed is formed, the melt temperature 1s 





t 
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decreased continuously for a period of 48 hours. At the end of this period 
the crystal on the probe should occupy about one-third of the volume of the 
jar. During this same period the probe and air temperatures will also de- 
crease, but as long as the initial equilibrium temperatures have been 
established the crystal on the probe should grow satisfactorily. However, 
as the melt temperature is decreased, the jar walls become cool enough to 
allow crystals to form on the inside of the walls. 

The melt temperature is decreased for another 24 hours. During this 
period the probe and wall crystals will start to meld into one another. 


When this occurs, the melt heater control no longer has much effect on the 


_erystal growth and the crystals will grow without any further temperature 


changes. It takes approximately another 48 hours for crystallization to be 


| Complete throughout the jar. Hence it takes about five days to complete the 


crystal-growing run once the seed is formed. Obtaining a good seed crystal 
is the most critical part of the procedure. If more than two or three 
Single-crystal seeds are allowed to grow, the individual crystals obtained 
will generally not be large enough to use. 

The formation of crystals on the jar walls is undesirable. Many dif- 
ferent temperature combinations were tried to prevent this from happening. 
However, whenever the jar walls were heated enough to prevent crystilliza- 
tion, the crystal on the probe would be limited to a smaller size than 
obtained with wall crystallization. Thus, the wall crystallization had to 
be accepted as one of the limitations of the apparatus. In a few of the 
runs, the probe and wall crystals melded into one large single crystal which 


Occupied almost the entire volume of the jar, but generally the wall and 


probe crystals would be separate single crystals. Even then, it was not 
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unusual to obtain a large piece of single crystal from the wall crystalliza- 
ion. 

To prevent the wall effect the apparatus would have to be modified in 
two main ways. First, the volume of the jar should be reduced to about 
half its present size in order to make the melt temperature controller ef- 
fective throughout the growing process. Second, the jar walls should be 
thinner. The walls of the present apparatus are relatively thick (- > in.) 
which causes a great deal of heat to be dissipated in the walls -- at some 
point, enough to sustain crystallization. Thinner walls would be more 
responsive to the melt temperature control and would require less heat to 
warm them sufficiently to prevent crystallization; thus, there would be 
less effect on the probe crystal. These changes were not accomplished 
because the crystals which were grown were large enough for the purposes 
of this work. 

Once crystallization was complete throughout the jar, the entire jar 
is removed from the refrigerator and transferred to the glove box. The 
procedure used to separate the single crystals from the solid mass of ben- 


zene is described in Appendix C. 
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Appendix C 


Crystal Cutting Methods 


Two methods were used to separate: the single crystals from the frozen 
mass which resulted from the crystal-growing procedure described in Appen- 
dix B. First, it was found that if the frozen mass were eee inside the 
cold box (- 10°C) for about an hour, sublimation along the crystal boundaries 
would occur more rapidly than from the crystal surfaces. As a result the 
boundaries become clearly visible and if a thin stainless steel spatula is 
inserted into a boundary at two or three points, the entire piece of single 
erystal would break away from the mass. However, if the crystal does not 
break off easily, this procedure should not be used, because additional 
force could possibly produce strains in the crystals. 

The second and more common method was to use a hot-wire cutter to cut 
out the single crystal pieces. Previous investigators have used sharp 
knives and even hacksaws to cut organic crystals, but it was feared that 
such a procedure would lead to severe straining of the crystals. The hot- 
wire cutter was made in the form of the miniature jig saw shown in Fig. C-l. 
The base is made of plywood and covered with a piece of plate glass. A 
Michrome wire (number 18 gauge), held taut: in a teflon collar by a large 
cantilever leaf spring, passes through the base to an electrical connection 
on the back. The other end of the wire goes from the teflon collar, over 
the spring, and to another electrode on the back side. The large wire size 
is used to make the cut wide enough to preclude refreezing the two pieces 
' back together once the cut is made. The two electrodes are connected to a 
Variac and just enough current is passed through the wire to melt the crystal 


Slowly as the crystal is pushed through. The crystals were cut into cubes 
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with sides 3.5 to 4.0 cm in length. Smaller pieces or odd shaped ones can 
also be used if the initial orientation procedure described in Appendix D. 
reveals that the crystal is oriented close to the b-axis. 

To prepare a piece of crystal for optical crientation, two flat, 
parallel sides must be melted onto the sample. This is done with the ver- 
tical planing apparatus shown in Fig. C-2. The vertical planer is the same 
pieces of apparatus shown in Fig. 7a except that the sample holder housing 


has been removed and the track is held in the vertical position. 
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PEC HOt wire Cutter 





Fig. C-2. Vertical planing apparatus. 
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Appendix D 


Optical Orientation Procedure 


A Vickers' M-72C Research Polarizing Microscope was used in conjunction 
with two multiaxes universal stages to accomplish orientation of the crystal 
(see Section 3.4). Coarse orientation was performed with a locally manu- 
factured three-axis universal stage which was patterned after a similar 
device used by Heseltine and By eae and which is shown in Fig. D-l. The 
three~axis universal stage or the “rough universal stage" (R.U.S.) consists 
of two gimbal rings which are mounted inside a piece of large diameter 
aluminum tubing and which rotate about mutually perpendicular horizontal 
axes. Another ring is set inside the inner gimbal and is free to rotate 
about an axis perpendicular to the plane of that gimbal. Once coarse ori- 
entation is obtained, the cyrstal is mounted in a 5-axis Leitz Universal 
Stage (L.U.S.). Each gimbal on the L.U.S. is graduated in one-degree 
increments which allows fine alignment of the crystal to be accomplished. 
A 10x objective lens is used during the coarse orientation and a 30x ob- 
jective is used with the L.U.S. The R.U.S. is used because it can accom- 
modate a larger piece of crystal than the L.U.S. Once the coarse orientation 
procedure is complete, enough crystal has usually been melted away to allow 
it to fit into the L.U.S. The microscope with the L.U.S. attached is shown 
in Fig. D-2. (Note: It is convenient to designate the directions of the 
vibration planes of the polars when crossed at 90° as N-S and E-W. Thus, 
the crosshairs in the ocular and the horizontal axes of the universal stage 
can also be described as N-S and E-W.) 

In the procedures described below the terms isogyre, optic axis figure, 


and acute bisectrix figure are used frequently. Although briefly described 
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Fig. D-l. Three-axis universal stage. 





Fig. D-2. Vickers' M-72C research polarizing microscope 
with Leitz 5-axis universal stage. 
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in Section 3.4, a much more detailed description of these terms and the 
interference patterns associated with them is given by Wahlstrom”’ and 
one. These references should be studied prior to attempting the 
orientation procedure. 

The first step in the optical orientation procedure is to take a 
relatively large piece of single crystal (3.5-4 em cube) and melt two 
parallel faces on it with the vertical planing apparatus described in 
Appendix C (see Fig. C-2). After the surfaces are prepared they quickly 
become opaque because of sublimation. To keep the faces optically clear, 
glass microscopeslides -- initially at room temperature -- are frozen onto 
each face. These slides also provide a means for suspending the crystal in 
Ene center of the R.U.s. 

The R.U.S., holding the randomly oriented crystal, is put onto the 
microscope stage. (The microscope is located inside the refrigerated glove 
box at this point.) The gimbals of the R.U.S. are rotated until maximum 
extinction occurs. (There is a possibility of seeing an isogyre during 
this observation; however, usually extinction will occur.) Once the maximum 
extinction position is found, the extinction directions, which correspond 
to the polar axes, are marked on the glass slide. The R.U.S. is removed 
from the microscope stage and two new parallel surfaces are melted onto the 
crystal, parallel to one of the extinction directions. The orientation 
procedure is repeated. If an isogyre is not seen now, a third pair of 
parallel faces is melted onto the crystal, parallel to the other extinction 
direction previously determined. As discussed in Section 3.4, one of the 


optic axes should be observed during one or the other of these observations. 
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Once an optic axis is found, the crystal is transferred to the L.U.S. 
To ensure that the L.U.S. is centered with respect to the microscope stage, 
the L.U.S. is fastened onto the microscope stage with two set screws. The 
holes in the base of the L.U.S. and in the micrcscope stage are aligned so 
that the L.U.S. will be centered when the set screws are secured. 

The next eight steps in the orientation procedure will be explained 
with the aid of Fig. D-3. A schematic drawing of the position of the crystal 
relative to the L.U.S. is shown for each step. (Note the drawing is two- 
dimensional whereas the actual L.U.S. has three degrees of freedom.) The 
top and bottom faces of the crystal as they appear in step l have been 
designated A and B, respectively, and maintain the same identity through- 
out the procedure. A top view of the glass slide which is frozen onto 
side A and/or B is shown with the position of the reference marks used 
in each step. 

Two reference marks are used. One coincides with the zero index mark 
on the inner ring of the L.U.S. and is used to keep the crystal in the same 
relative position with respect to the ring throughout the procedure. This 
mark will be called the “index reference mark". (Note: The crystal can 
be held secure in the L.U.S. by first putting a thin coat of vacuum grease 
on the inner ring and then pressing the glass slide, which has previously 
been frozen onto the sample, into place.) The second reference mark is 
made in a N-S direction and is used to ensure that the relationship between 
the N-S and E-W directions of the crystal with respect to the horizontal 
plane of the L.U.S. in step 8 is the same as in step l. If it is not, one 
can easily become confused and there is a possibility that subsequent cor- 


rections will be made in the wrong direction. 
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The procedure is as follows (see Fig. D-3): 

Step 1 -- The crystal is suspended in the L.U.S. by the microscope 
slide attached to side A and the combination is put onto the microscope 
stage. The gimbals of the L.U.S. are rotated to bring one of the optic 
axes as close to the center of the field of view as possible and to make 
the isogyre parallel to one of the crosshairs in the ocular. The N-S and 
the index reference marks are put on slide A as indicated. (The initial 
direction of the N-S reference mark is arbitrary but once a direction is 
chosen, it must be maintained throughout the sequences. The N-S index mark 
should not be made until the desired pattern is centered in the field if 
view). 

Step 2 -- The L.U.S. is removed from the microscope stage and the 
crystal is set on top of the L.U.S. in the same relative position as it 
had while suspended. The N-S reference mark is unchanged. The index ref- 
erence mark is transferred to slide B. 

Step 3 -- The N-S reference mark is transferred from slide A to 
slide B. Slide A is removed by pressing a small flat piece of metal, 
which is at room temperature, against the slide; the heat transfer will melt 
the bond between the slide and the sample in about 2-3 seconds. 

Step 4 -- A horizontal surface is melted onto side A with the vertical 
planer; then a glass slide is attached to the new surface and the N-S ref- 
erence mark is transferred from slide B to slide A. 

Step 5 -- The N-S and E-W gimbals on the L.U.S. are set to their zero 
positions. Reference marks are unchanged. 

Step 6 -- The crystal is turned over so that side A rests on top of 


the L.U.S. The same N-S reference mark is on slide A as in steps 4 and 5; 
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ILLUSTRATION OF Steps Usep DURING ORIENTATION PROCEDURES 


Front View of Crystal on Top View of Reference Marks on 
Universal Stage (U-stage) Glass Slides (A or B as indicated) 


STEP 1 


A 


Index 


E N E 
B 


Slide A 





US. 
Crystal is suspended in L.U.S. in an orientation at which on isogyre, optic 
axis pattern, or an acute bisectrix pattern has been observed with the 
microscope. If orientation is close to b axis, angles have been corrected 
for refraction. Reference marks on slide A in the N direction and in line 


with index on inner ring of the L.U.S. 


STEP 2 


Crystal is set on top of L.U.S. in same relative position as it had when 
suspended. Reference mark on slide A in N direction. Index reference mark 


s transferred to slide B. 


Figure D-5 
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STEP 5 





A 
ty 


N-S reference mark from slide A is transferred to slide B. Index reference 


mark unchanged. Slide is taken off side A. 


Step Y 
ES 


A new horizontal surface is melted onto side A with the vertical planer. 


H 
w > 


A slide is frozen onto the new surface and the N-S reference mark is trans- 


ferred to slide A. 
A STEP 5 
ANTA 
B | 


The N-S and E-W gimbals on the L.U.S. are set to zero. Reference marks 


are unchanged. 


Ficure D-3 (cont’D), 








Crystal is turned over so that side A rests on top of the L.U.S. The same 
N-S reference mark is on slide A as in steps 4 and 5; however, now arrow is 
pointing in the S direction. Index reference mark is transferred from 


Serlesprto slide A. 


STEP / 


Bez 


Denet aken off of side B and a horizontal face is melted onto the 
crystal with the vertical planer. Another slide is then frozen onto new 


surface. Reference marks are unchanged. 


STEP 8 





Crystal is turned over and suspended in the L.U.S. Side A is on top and in 
same relative position as in step 1. Crystal is ready for another obser- 


vation with the microscope. 


Fieure D-3 (CONT D): 
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however, now the arrow is pointing in the opposite direction. The index 
reference mark is transferred from slide B to slide A. 

Step 7 -- The slide istaken off of side B and a horizontal face is 
melted onto the crystal with the vertical planer. Another slide is frozen 
onto the new surface on side B. Reference marks are unchanged. 

step 8 -- The crystal is turned over and suspended in the L.U.S. 

Side A is on top, in the same relative position as in step 1. The crystal 
is now ready for another observation with the microscope. 

During the next observation, the optic axis figure should be close to 
the center of the field of view with the N-S and E-W gimbals positioned at 
lesser angles than those observed before the correction was applied; if a 
perfect correction was made, the N-S and E-W gimbals would be in their zero 
positions. If the optic axes is still far off the zero position (greater 
than 10°), it should be centered once again and steps 1 through 8 should be 
repeated. 

When an optic axis can be seen close to the zero position, the inner 
ring of the L.U.S. should be rotated until the isogyre is parallel to one 
of the crosshairs in the ocular. Next, the N-S or E-W gimbal (the one used 
depends on the direction of the isogyre) is rotated so as to move the optic 
axis along the crosshair in the ocular which is parallel to the direction 
of the isogyre. As the optic axis is moved, the field of view will either 
become dark and the isogyre will broaden EE Or the Acute bisectrix 
figure and the other optic axis figure will be seen in sequence. (See 
Fig. 5 in Section 3.4,) If the former is observed, the gimbal must be 
rotated in the opposite direction to obtain the acute bisectrix figure. 


Once it is observed, the gimbals should be rotated until the cross coincides 
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with the crosshairs in the ocular. In this position the acute bisectrix, 
1.e., the D axis, is aligned parallel to the vertical axis of the micro- 
scope, but the N-S and E-W gimbals are at some angle with respect to the 
horizontal plane of the microscope; hence, the observed angles must be cor-- 
mectea tor refraction, 1.e€., Snell's Law, to obtain the true angular rela- 
tionship of the crystal with respect to the horizontal plane. Since the 
light is propagating close to the b-axis at this point, the index of refrac- 
tion (1.64) corresponding to that direction is used. Y (The observed angle 
must always be reduced to obtain the true angle.) After the corrections 
have been made, steps 1 through 8 are repeated except that now, in step 1, 
the acute bisectrix figure has been centered (and corrected for refraction) 
rather than the optic axis figure. The procedure of observing, correcting, 
and preparing new surfaces (steps 1-8) is repeated as often as is necessary 
to bring the b-axis coincident with the vertical axis of the microscope with 
the N-S and E-W gimbals in their zero positions. If the N-S reference marks 
are in the same direction after each sequence, each succeeding observation 
will require rotation of the N-S and E-W gimbals in the same direction as 

in the previous observation but the amount of correction required will be 
less each time. If the angles become greater, the "corrections" have been 
applied in the wrong direction. If one carries the procedure out carefully, 
an orientation to within +2° can be obtained. 

A practical difficulty in these steps arose because the cold box was 
not large enough to accommodate both the microscope and the planer surface 
Simultaneously; hence, each SÉ of apparatus had to be alternately taken 
out and put into the box to accomplish the required EE A larger 


work area would have been desirable. 
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If the b-axis is the axis to be oriented, the orientation is maintained 
with the glass slide which has been frozen onto the final surface. To 
orient the a or c axis the appropriate isogyre is aligned with the 
E-W crosshair in the ocular. The crystal is then set on top of the inner- 
most ring and the index reference mark, which was previously put onto the 
upper slide, is transferred to the lower slide with the two marks agreeing 
as closely as possible. The L.U.S. is returned to the microscope stage and 
the alignment is checked with the 3x objective lens. The pattern seen with 
the 3x objective is much more diffuse than the corresponding pattern seen 
with the 30x objective; consequently, it takes some practice to correlate 
the two patterns correctly. By alternately viewing the crystal in the sus- 
pended position with the 30x objective and in the "on top" position with 
the 3x objective, it was found that if the procedure is carried out care- 
Fully, a correction of a degree or two in the horizontal plane was all that 
was necessary to compensate for the movement. 

Once the crystal is oriented in the “on top" position, the L.U.S. is 
transferred to the planing apparatus shown in Fig. D-4. The L.U.S. is 
fastened onto the platform of the planer in the same relative position as 
it had when it was on the microscope stage. Thus, the E-W position on the 
platform coincides with the axis aligned on the microscope stage and the 
face of the heated planer is exactly perpendicular to the axis. Consequently, 
a face can be melted onto the crystal which is perpendicular to the desired 
axis. The final surface is sealed with a glass slide. 

The orientation of the a or c axis can be checked again with the 
L.U.S. and the 30x Se lens. The directions of the a and c axes 


correspond to the obtuse bisectrix and to the optic normal of the indicatrix, 
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respectively. Each of these directions also has a characteristic inter- 
ference pattern associated with it (described by Wahlstrom” 9, However, the 
resolving power of the 30x objective lens is not great enough to obtain as 
well defined patterns as seen along the acute bisectrix. Consequently, 
these patterns were used only as a secondary check on the alignment. The 
resolving power of the 30x objective also limited the amount of the acute 
bisectrix pattern which could be seen in the field of view. Usually the 
optic axes would be at the extreme edge of, or just outside, the field of 
view. Therefore, to determine which isogyre corresponds to the a axis, 
the crystal has to be rotated back and forth on the N-S and E-W gimbals to 
see the optic axes. After the a axis has been identifed, the cross can 
be centered in the field. It would be desirable to use a 40x or 50x objec- 
tive lens with the L.U.S., to obtain a greater resolving power which, in 


turn, would yield sharper interference patterns. 





Fig. D-4. Planing apparatus used for preparing surfaces 
perpendicular to the a or c axis. 
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Appendix E 


Sample Cutting and Polishing Procedure 


To cut a sample on the apparatus described in -Sec. 3.5, the 
apparatus is used in the vertical position as shown in Fig. 6b, The cavity 
is removed from the sample holder housing and a small glass plate ( 3 x 3 
inches) is set on top of the upper face of the housing. A piece of oriented 
crystal is frozen onto the glass plate with a few drops of benzene, using 
the microscope slide which was attached during the orientation procedure 
as a base. Hence, the oriented axis is perpendicular to the glass plate 
and parallel to the axis of the cutter. The sample is cut by bringing the 
heated cutter down onto the sample for 4-5 second intervals. The cutter 
cannot be kept on the sample continuously because the sample will be melted 
into the shape of a truncated cone by the heat inside the cutting cylinder. 
Therefore, the cutter is raised off the sample after each cutting interval 
to ensure that the sample will be cylindrical and also to reheat the tip 
of the cutting tool. Thus, the cutter is alternately raised and lowered 
until a sample is cut out. It was found that the use of two cutting tools 
with diameters of 7/8-inch and 3/4-inch, respectively, produced a more 
cylindrical sample than obtained with making only one cut with a 3/4-inch 
diameter tool. After each cut is complete, the excess crystal is trimmed 
away with a sharp razor blade. The result is a sample which is approxi- 
mately 3/4-in. (vw 1.9 cm) in diameter and one inch long; the glass slide 
is still frozen onto one face of the sample to preserve orientation. 

Nica ample on the microscope slide, a brass spacer plate 
EE Elek Beste thick, and has a 3/4-in. diameter hole in 


the center is slipped over the sample so that one face of the plate rests 
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on the slide. The sample holder is set on top of the spacer plate and the 
sample is frozen to the holder by inserting a few drops of liquid benzene 
into the small space between the sample and the holder. The oriented 
crystal axis is now parallel to the sample holder axis. In order to be 
able to prepare flat, parallel faces on the sample, approximately 1/8-in. 
of benzene must protrude beyond each side of the sample holder. The spacer 
plate provides the necessary excess on one side whereas the amount of excess 
on the opposite side is determined by the length of the piece of crystal 
being prepared. Thus, the length of the oriented piece of crystal must be 
at least 1/4-in. greater than the thickness of the sample GEES 

The sample holder and spacer plate along with the enclosed sample and 
the attached microscope slide are mounted on the planing apparatus shown in 
Fig. 6a. The holder is held flush against the cavity in the sample holder 
housing by two long screws which pass through the housing and the cavity and 
thenscrew into the holder. Hence, the face of the planer and the side of 
the sample holder are exactly parallel. The glass slide and the spacer 
plate are removed and the heated planer is used to melt a flat face onto 
the sample. The sample is then turned around, secured once again to the 
cavity, and a face is melted onto the opposite side of the sample. The 
temperature of the planer surface is maintained about 2~-3°C above the 
freezing point of benzene during the facing operations. The sample faces 
are melted to the point where approximately 1/16 of an inch of benzene still 
protrudes beyond the sides of the sample holder. The excess benzene is 
rubbed off during the final polishing procedure. 

Many different types of abrasives and polishing cloths were tried in 


efforts to make the final surfaces of the sample smooth and flat. Because 
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benzene is relatively soft, it was found that even the finest grit abrasive 
Papers left scratches on the crystal surfaces. A velveteen polishing cloth 
was finally used to obtain the requisite surfaces. (The velveteen was 
procured from the J. J. Morris Company, Southbridge, Mass.) The velveteen 
cloth has an adhesive backing; consequently, a polishing ‚surface was made 

by attaching it to a flat glass plate. To obtain the final surfaces on the 
crystal, the sample holder is set onto the plate and the excess benzene 
rubbed off by pushing the sample holder across the velveteen surface in one 
direction with short, positive strokes. At the end of every third or fourth 
stroke, the sample holder is turned 90°. This procedure is repeated until 
the benzene is flush with the flat side of the sample holder. The paral- 
lelism of the sample is determined by the degree of parallelism to which 

the sample holder could be machined. For this case, the sides of the sample 
holder are parallel to within 0.0005 in. 

After the final surfaces of the sample have been prepared, one face is 
covered with a piece of 1/4-mil clear mylar which has been cut into a 1.25-in. 
diameter disc. The edges of the disc are sealed onto the sample holder with 
a few drops of a glycerine-water solution (to be described later); then the 
mylar is covered with a teflon plate which is screwed onto the sample holder. 
The mylar-teflon seal prevents sublimation Serstbesy ep. Je he 
opposite sample face is covered with a transducer mylar seal which is 


described in Appendix F. 
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Appendix F 


Bonding Procedure 


The bonding procedure is initiated by cutting a piece of 1/4-mil 
aluminum coated mylar into a 1.25-in.diameter disc and removing the oxide 
layer from the aluminum by immersing the disc in concentrated nitric acid 
and then rinsing with distilled water. The transducer is bonded onto the 
mylar at room temperature with Nonaq stopcock grease. Prior to making the 
bond, a small amount of Nonaq is heated on a glass slide to boil off any 
absorbed water vapor. After the Nonaq is allowed to cool, a drop is put 
onto the mylar and the transducer is bonded by twisting it back and forth 
on the mylar with considerable thumb pressure until all the excess grease 
is squeezed out of the bond. It is important to note that the transducer- 
mylar bond must be made at room temperature in order to obtain a thin, strong 
bond. If the procedure is carried out at - 10°C, a poor bond will result, 
especially when utilizing thick transducers, (During the course of the 
experiment, a peculiar property of the Nonaq stopcock grease was observed. 
It was found that Nonaq which had been allowed to set and absorb water vapor 
from the air, and subsequently heated to boil off the water vapor, yielded 
a much stronger bond than Nonag from a freshly opened bottle. The "old" 
Nonaq had a milky, opaque appearance and a smooth consistency whereas the 
"new" Nonaq appeared to be translucent and had a grainy consistency. There 
was a considerable difference in the bonding properties of the two types.) 

Once the transducer-mylar bond was made, the combination is cooled down 
in the cold box and,immediately after one face of the prepared sample is 


sealed with the teflon-mylar combination, the transducer-mylar combination 





Zuge 


is bonded onto the opposite face with a glycerine-water solution. It was 
found that if pure glycerine was used, the glycerine would freeze onto the 
benzene before a thin bond could be made, which resulted in a very poor 
ultrasonic echo pattern. After a great deal of experimentation, it was 
determined that a solution of 90% glycerine and 10% distilled water by 
weight exhibited the best bonding characteristics. The bond would remain 
fluid at - 10°C and would freeze at approximately - 20°C, 

To bond the transducer-mylar combination onto the benzene, a small 
drop of glycerine solution is put onto the center of the sample face and 
the transducer-mylar combination is set on top. The ends of two Q-tip sticks 
are used to twist the transducer-mylar combination back and forth gently 
until the glycerine is spread evenly across the sample face and the excess 
is squeezed out. The aluminized mylar seals the sample and provides the 
ground connection for the transducer. 

The experimental limitations of the bond are discussed in Section 3.6. 


An estimate of bond losses and how they were obtained is given in Section 4.2, 
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Appendix G 
Apparent Sound Absorption Due to Transducer Thickness 


The results of oy ese mame” experiments in determining the sound 
absorption due to differences in transducer thicknesses were briefly pre- 
sented in Section 4.2. Although Sylwestrowicz offered no explanation for 
the thickness effect, one possible cause of the effect could be the change 
in the quality factor (Q) of a transducer when it is operated above its 
fundamental frequency. The Q of a transducer depends on the specific 
acoustic impedences of the transducer material and that of the substance 
on which the transducer is bonded. Foran air-backed transducer the Q is 


75 
defined by 


S T 
Q zn 5 Go (G.1) 


Here n is the number of the harmonic at which the transducer is being 
excited and (ps), and (ps), are the specific acoustic impedences of the 
transducer and of the solid, respectively. Equation G.1 shows that the Q 
of the transducer depends on the loading conditions imposed by the sample. 
In general, the presence of the bonding agent will cause the Q to be 
lower than predicted by Eq. G.1. 

For a given loading condition the Q of the transducer can be defined 
in terms of its bond width characteristic as 


a 
Q = Av S (G.2) 


Av is the bandwidth of the device at the resonant frequency Ya and is 
defined by the points at which the acoustic power is 3 db (50%) below the 


peak power of the acoustic pulse at the fundamental resonant frequency. 





a= 


o 


Hence, 1f = = 10% at the fundamental frequency, the bandwidth is 500 KHz 
for a 5 MHz transucer whereas it is 1500 KHz for a 15 MHz transducer. The 
value of Av is a characteristic of a given device; consequently, does not 
change with the harmonic at which the transducer is operated. Therefore, 
for a 5 MHz transducer operating at 15 MHz the Q is three times that of 
a 15 MHz transducer operating at its Fundamental, * 
A relationship between Q and the energy transmitted into the detec- 
tion system can be obtained by first considering the spectral distribution 
of the voltage response V(v) of the transducer. For a voltage which con- 
tains components in the frequency range Gë + v , the total voltage seen by 
the detection system is proportional to the area defined by the curve which 
describes the spectral distribution on a voltage vs. frequency plot and the 
two vertical lines drawn through the points tv on the frequency axis. If 


a Gaussian curve is assumed, the area can be expressed in terms of the 


Gaussian distribution function, which can be written in the form 


A vs? 
Area = Ee | e 2 S dv 
Yon 
V 
r 
1/Q 
2V 2 
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Since the error function (erf) decreases with decreasing values of its 
argument, Eq. 6.3 shows that the area is smaller for the system with the 


higher Q. The transmitted energy is proportional to the square of the 


ge 


voltage--hence, proportional to (Area)?. Hence, for a pulse of acoustic 
power which has a given frequency spectrum, the device with the lower Q 
will pass more of the frequency components into the detection system which, 
in turn, yields a lower apparent attenuation. 

If Eq. G.3 is used to calculate the area for two different values of 
Q that have a ratio of 3:1, it is found that the system with the higher Q 
transmits 81.5% less energy than the system with the lower Q. Sylwestrowicz 
made sound absorption measurements on three different materials at 15 MHz 
with transducers which had fundamental frequencies of 5 MHz and 15 MHz. 

The apparent attenuation obtained with the 5 MHz transducer was 73%, 88%, 
and 25% greater for the three substances, respectively, than that obtained 
with the 15 MHz transducer. Hence, there is a qualitative agreement between 
the predicted and observed effect of Q on the apparent attenuation in two 
of the three cases. 

The procedure described by Sylwestrowicz was used to obtain an esti- 
mate of the apparent attenuation due to the thickness effect for this experi- 
ment. Transducers with fundamental frequencies of 2.0, 2.5, and 5.0 MHz 
were successively bonded onto the same sample (oriented along the c-axis) 
and absorption measurements were taken at 250°K. The results obtained with 
each transducer is shown in Fig. G-l. The several measurements exhibit the 
same frequency dependence but are displaced with respect to each other in 
magnitude. Next, the absorption values at a particular frequency, e.g., 5, 
6, 7.5 MHz etc., are taken from each curve and are replotted as a function 
of transducer thickness. This generates another family of curves (Fig. G-2) 
on an absorption versus transducer=thickness plot,which are parametric in 


frequency. It is assumed that the absorption associated with transducer 
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thickness vanishes at zero thickness. Hence, the curves of absorption versus 
transducer thickness are extrapolated to zero thickness and the difference 
between a at the intercept and the value of a particular transducer thick- 
ness represents the energy loss related to the transducer at the frequencies. 
of interest. By applying the corrections obtained from Fig. G-2 to each 
experimental point obtained with the three transducers, the three curves 
collapse into the single corrected curve shown in Fig. G-1. A comparison 
of the measured and corrected absorption values is made in Table G-1. It 
can be seen that the corrections for the 5 MHz transducer are much less 
than for the thicker 2.0 and 2.5 MHz transducers. 

The experiment did show the existence of an apparent attenuation due 
to transducer thickness; however, further experimentation would be necessary 
to be able to make corrections to all the data obtained. Transducers with 
10 and 20 MHz fundamentals should be used to obtain more points on the ab- 
sorption versus transducer-thickness plot which, in turn, would make the 
extrapolation to zero thickness more accurate. Also, values should be ob- 
tained as a function of temperature for each transducer to determine the 
temperature dependence of the effect. Such an experiment would be a for- 
midable task and was not carried out in this work due to the difficulties 
associated with the sample preparation. It would be feasible, however, to 
perform the experiment with a "hard" solid with moderate intrinsic attenua- 


tion. 
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Table G-1 
Correction of a Due to Transducer Thickness Effect 


eta ola Ponto cto ae ale 
dp dp dh dp 


v E correction ` ` Baken % 
(MHz) (cm!) (cm 1) (em!) correction 
5 0.156 07023 0.133 14.7 
6 02235 0.072 0.160 31.8 
Veo 0.294 0.081 0.213 21.9 
10 0.433 0.118 OSS 2722 
Do 0:3923 0123 0.800 2325 
14 0.807 0.367 0.440 OO 
15 0.604 0.124 0.480 20% 3 
17,5 0.905 0.305 0.600 33:8 
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